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ERRATUM 


Two errors the Appendix have been pointed out the author the paper 
radio waves ionized gas thermal equilibrium” (Can. Phys. 38, 1114, 1960). These 
errors not occur the corresponding analysis Pines and Bohm, which correct within 
its range validity. 

The first error consists the omission the second order term the expansion the 
logarithm leading equation (58). The second occurs the final integration leading 
equation (73), where the lack independence Fourier components with wave number 


vectors equal magnitude and opposite direction was not taken into account. 

The corrected forms equations (6), (7), (73), and (76) are obtained replacing 
the original equations The corrected form equation (8) obtained replacing 
the original equation The results obtained for the frequency spectrum scattered 
power are not affected. 
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THE CRACKING ACTIVITY ICE DURING CREEP! 


ABSTRACT 


Observations were made one temperature the cracks that form ice 
during creep under constant compressive load. The ice had hexagonal symmetry 
with respect the grain boundaries. The load was applied perpendicular the 
long axis the grains. piezoelectric crystal was used detect the cracking 


activity the ice. 

Two stages cracking were observed. The first occurred during the transient 
period the creep, and the plane these cracks tends parallel the grain 
boundaries and the direction the stress. The rate which these cracks 
formed decreased very markedly the creep rate approached constant value. 
Above certain stress, the creep rate continuously increased with time. Under this 
condition the second stage cracking was observed. These cracks tend 
more irregular direction and occur planes that are 45° the applied 


stress. 


INTRODUCTION 


Observations the elastic and creep behavior ice under constant com- 
pressive stress have been carried out Brown (1926). lower stresses 
observed creep rate which, after the initial transient stage, tended decrease 
with time constant value characteristic the stress. critical stress 
range, between 190 and 200 p.s.i. for the temperature and load conditions 
the experiment, the character the curve changed radically from 
one for which the creep rate decreased with time one for which increased 
with time. The period observation was never more than hours. 

Throughout his paper, Brown refers internal cracking the ice and 
associated audible that occurred during load test. Qualitatively, 
his comments indicate that the intensity depends the tempera- 
ture and the magnitude the applied stress. 

The creep ice and its general dependence stress, temperature, and time 
similar that observed for many materials. Glen (1958), his review 
such investigations, found that the minimum observed strain rate related 
the applied stress according the law: 


(1) 
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where the strain rate, the applied stress, and constant. The 
exponent observed increase with within the range 1.85 4.16 for 
1.0 15kg/cm?. Glen (1955) observed that until accelerating creep occurs, 
ice obeys the Andrade law with the initial transient the creep curve being 
proportional the cube root the time. 

Cracking activity has been reported other materials during compressive 
loading. Obert (1945) conducted investigations the laboratory see the 
microseismic activity, which observed rock the field, could also 
detected with rock specimens subject compressive loads the laboratory. 
observed high ‘‘microseismic when the load was between 1/8 1/2 
the ultimate strength the rock. This rate had decreased minimum when 
the load reached 1/2 3/4 the ultimate strength. the load was increased, 
observed that the rocks that tended shatter showed marked increase 
activity about 8/10 the crushing strength, whereas 
the rocks that tended crush showed comparatively little increase activity. 
was able show that the first maximum his observed 
activity was due the geometry loading. held the load constant 
stresses 9/10 the ultimate load, the activity rate the sample decreased 
with time constant value which depended the magnitude the load. 

Bridgman (1949) attributed anomalies the volume change during loading 
for loads approaching the ultimate strength the material the forming 
internal cracks interstices grain boundaries, mechanical impurities, 
around dislocations the lattice. 

Jones (1952) observed that the velocity ultrasonic pulse transmitted 
through concrete specimen direction perpendicular the applied com- 
pressive load began decrease when the load was about 30% the 
ultimate strength. attributes this decrease the formation internal 
cracks. Because the velocity pulse propagated the direction the stress 
did not change, concluded that the plane the cracks must parallel 
the direction the stress. 

Riisch (1959) using crystal pickup, recorded the internal cracking 
concrete subject compressive loads. For loads within 85% the 1-minute 
ultimate strength, observed high rate cracking when the load was first 
applied. This rate dropped off but increased again rapidly just before failure 
occurred. 

The purpose the observations reported this paper was investigate 
one temperature and under constant compressive load, the cracking activity 
that occurs ice during creep. particular, was desired investigate 
the sudden transition that Brown reported the character the 
curve and see this transition was some way related the internal 
cracking. 


EXPERIMENTAL PROCEDURE 


The observations were carried out rectangular ice specimens 5X10 
and long prepared —10° with standard metal cutting lathe and 
wood planer. small, manually driven testing machine (Hounsfield 


an — 


PLATE 


Fic. Thin sections test specimen viewed with polarized light: section perpendicular 
direction growth, section parallel direction growth. 
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tensometer) with capacity 2000 was used for applying the load. The 
test pieces were mounted between special loading heads one which shown 
Fig. The silicone grease, separated from the ice the silicone rubber 


DIRECTION APPLIED STRESS 


ORIENTATION GRAINS 


TEST PIECE 
LONG 


SILICONE RUBBER 
OIAPHRAGM SILICONE GREASE ———~ 


Fic. Sketch loading head construction. 


diaphragm, ensured that the applied compressive stress was uniform over the 
test piece. 

keep the silicone rubber from shearing during experiment, water was 
frozen the small space between the test piece and the loading heads. 
The frozen bond between the ice and failed soon load was applied 
but the ice was constrained from deforming laterally. This constraint influenced 
the distribution cracks the vicinity the loading head will shown 
later. 

extensometer with 2-in. gauge length and strain magnification 3300 
was used for measuring the creep. The gauge had remounted periodically 
because its limited range. The creep was recorded Baldwin X-Y 
recorder, model The pen the recorder automatically and continuously 
traversed the chart paper constant rate inches per minute. The 
curve was constructed aligning the measured sections the 
known time base. 

The ice used the observations was prepared from ordinary tap water 
cylindrical container about meter diameter. The direction freezing was 
perpendicular the water surface. When the water surface cooled ice 
grains were sprinkled the surface act centers nucleation. The result- 
ing grain structure seen with polarized light shown Fig. Figure 
thin section cut parallel the ice surface and Fig. section cut perpen- 
dicular the surface. The columnar structure typical for multigrained ice 
often encountered under natural conditions. 
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The test specimens were cut that the long axis each grain was perpen- 
dicular the face and therefore perpendicular the applied 
compressive stress shown Fig. The stress applied the test specimen 
was kept constant +1% the applied load during the test. 

measure the cracking activity the ice during creep, piezoelectric 
crystal was frozen one face the test piece. The output from the crystal 
was fed into high-gain amplifier and then recorded. Each time disturbance 
occurred within the ice that would set shock wave, the event caused 
deflection the recorder pen. 

The observations were carried out cold room maintained 
The test pieces were stored the room kerosene for least hours before 


they were loaded the testing machine. 


OBSERVATIONS 
Cracking Activity 

Figure shows that the observed curves for all tests exhibit 
initial increase creep rate with time. Those subjected stress less than 
about had decreasing creep rate with time after about the first 
half hour. Two these showed tendency later period increasing 
creep rate with time. 

Figure shows the observed dependence the total cracking activity 
time for three test pieces subject the applied stress given the figure. Also 
shown are the associated curves. stress kg/cm?, very 
little internal cracking stress 15.2 kg/cm? the cracking activity 
was appreciable. Most the activity occurred the first hour during the 
transient portion the creep-time curve. After hours, practically 
further cracking occurred. When the stress was greater than about 
and the creep rate showed tendency decrease with time, the cracking 
activity greatly increased. 

Some the details the curves are shown Fig. where 
the number disturbances that occurred 5-minute interval are plotted 
against time. The resulting curve approximation the rate cracking 
activity. One further set observations made stress 15.4 kg/cm? 
shown illustrate the oscillatory character the rate cracking activity 
after the first maximum. The cracks seem occur bursts though 
critical situation developed within the test piece, which then was relieved 
the formation cracks. Figure shows also that even when the applied stress 
was great enough that decrease was observed the rate creep, decrease 
was observed the cracking rate but then the rate increased beyond the 
original maximum. 

Fig. the logarithm the total cracking activity that occurred during 
the first hour, which amounted least 60% the total cracking activity 
observed, plotted against the applied stress. The line drawn the least 
squares fit all the observations made. This figure shows the marked depen- 
dence that the cracking activity has the applied stress. 


UNIT CREEP 


UNIT CREEP 


200 


100 


50 


200 


50 


Crack Formation 

The cracks were observed propagate parallel the long direction the 
grain boundaries and perpendicular the applied stress. This shown 
Fig. where each light mark crack. The edge view shows that some 
cases the cracks will propagate from face face. For stresses below about 
the crack distribution, which resulted during the period loading, 
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Creep-time curves. The stress kg/cm? given for each curve. 
Fic. Three typical examples the dependence the total cracking activity and the 
creep time the applied stress given. kg/cm?; (2) stress 
15.2 kg/cm?; (3) stress 16.1 
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Fic. Time dependence the rate cracking activity. 
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Fic. Stress dependence total cracking activity which occurred during the first 
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Fic. Examples crack development ice during creep (scale inches): stress, 
duration load, 180 minutes; stress, 15.2 duration load, 300 minutes; 
(c) stress, kg/cm?; duration load, minutes; (d) stress, 15.6 kg/cm?; duration load, 
225 minutes. 
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was fairly uniform shown Figs. and For stresses greater than about 
15.5 marked change occurred the character the cracking and 
this change was associated with the increase the cracking rate shown 
Fig. for specimen these cases, after initial cracking period 
minutes, during which the cracks formed uniformly throughout the test piece, 
the cracks tended form planes that were about 45° the applied stress 
and parallel the grain boundaries shown Fig. 7c. only one case out 
the five which these cracking planes developed was the plane not parallel 
the grain boundaries. This shown Fig. 7d; its curve 

Most the test pieces were quite free cracks triangular zone adjacent 
the loading heads shown Fig. 7c. This thought due the con- 
fining influence the loading heads. 

Two test pieces were subjected stresses the order kg/cm? and 
allowed deform until number cracks had developed. The angle between 
the plane each crack and the direction the stress was measured. The 
resulting distribution diagram for cracks observed, which shown 
Fig. shows that there definite preference for crack form with its 


NUMBER CRACKS 


ANGLE (DEGREES) 


Fic. Number cracks 2}-degree interval with plane crack angle stress. 


plane parallel the direction the applied stress. cracks were observed 
with their plane angle greater than 45° the stress. already mentioned, 
the long direction the cracks parallel the grain boundaries. 

Out cracks studied, occurred grain boundaries. additional 
eight occurred grain boundaries but were also partly transcrystalline. The 
remainder were completely transcrystalline. 

The water molecules ice are loose-packed and the single crystal has 
hexagonal symmetry. attempt was made see there was any correlation 
between the direction the axis hexagonal symmetry and the direction 


20 
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the plane those cracks that were transcrystalline. etching technique 
described Higuchi (1958) was used determine the crystallographic 
orientation. positive correlation was found because proved almost 
impossible form etch pit grains containing crack using Higuchi’s 
technique. 

was realized during the etch-pit observations and later confirmed 
observations using polarized light that there was bias the crystallographic 
orientation the grains. had been hoped that seeding with snow would 
result random orientation. Hillig (1958) has shown that ice crystals grow 
more quickly perpendicular the axis the hexagonal symmetry than 
parallel it. result, crystals with their symmetry axis perpendicular 
the direction growth tend squeeze out those that are parallel. the 
present experiments, crystals were observed that had their axis hexagonal 
symmetry within 45° the direction growth. Thus, the easy slip plane 
the grains tended parallel the grain boundaries. Because the symmetry 
axis tended normal the grain boundaries and therefore parallel the 
ice surface, was impossible with the apparatus available determine its 
direction accurately with polarized light. was observed that least 
the transcrystalline cracks were either parallel perpendicular the 
symmetry axis and therefore parallel perpendicular the easy slip plane. 


Deformation the Test Pieces 

Measurements were made one specimen that had been subject stress 
kg/cm? determine the change volume that occurred over the section 
where the creep was measured. increase volume approximately 
was observed. This increase probably due the formation the cracks. 
When the sample was immersed kerosene for storage after the observations, 
air bubbles were observed issue from some the cracks. 

All the lateral deformation the specimen was observed normal 
the grain boundaries. significant creep the direction parallel the long 
axis the grains was observed any the tests. The ice deformed way 
similar pile logs subject load perpendicular the long dimension 
the logs. the completion test, the edges most the test pieces 
were very irregular and was possible recognize individual grains because 
their displacement relative their neighbors. 


DISCUSSION 


earlier paper, Gold (1958) presented evidence which indicated that 
when ice, temperature between and —40° subject uniaxial 
compression normal the grain boundaries, there relaxation the shear 
stress the grain boundaries. Consider the array hexagonal grains shown 
Fig. the shear stress relaxes the grain boundaries, tensile stress 
will developed across the planes nearly parallel the applied stress, such 
whose average value given approximately 
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STRESS 


Fic. Array hexagonal grains. 


the average tensile stress, 

the compressive stress, 

proportion shear stress still effective the grain 
boundary, 


the angle between the applied stress and the plane under con- 
sideration. For 0°, 


This demonstrates that under conditions uniaxial compression, the 
shear stress can relax the grain boundaries, each grain will subject 
both compressive and tensile stresses. The magnitude these stresses 
each point will depend the value the applied stress, the local geometry 
the grain boundaries, and their orientation with respect the applied 
stress. the present case, the stress field biaxial because the columnar 
nature the grains. Under these conditions, one would expect cracks form 
parallel the grain boundaries with their plane tending parallel the 
direction the applied stress, that is, form planes across which tensile 
stresses can develop. Furthermore, would expected that the number 
such sites available for crack formation would limited and would 
depend the magnitude the stress was observed (Fig. 4). 

Figures and show that the initiation and extension crack ice 
time dependent. one considers the conditions necessary for the development 
crack according the theory put forward Griffith (1924), this would 
imply that (a) the stress across the rupture plane increases with time, 
imperfections the single crystal the grain boundaries coalesce until 
time nucleus the right size has formed from which crack can develop 
(Cottrell 1959), (c) both processes occur together. 

Cottrell states that for holes incoherent grain boundaries grow, the 
following relation must satisfied: 


(3) 


1146 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


where the tensile stress normal the boundary, 
the surface energy ice per unit area, 
the radius the hole. 


Figure shows that for the duration load associated with the observations, 
the compressive stress has greater than about kg/cm? for crack 
develop. Assuming that the shear stress has completely relaxed the grain 
boundary, this would give for value about kg/cm?. The value 120 
ergs/cm? for (Skapski al. 1957) gives, for the critical radius, 
cm. This value for little smaller than that used Cottrell calculations 
the critical tension for metals. 

Figure shows that the transient portion the creep-time curve occurs 
during the period when the cracking maximum. Glen (1955) has attributed 
this transient stage, and the condition increasing creep rate with time, 
recrystallization. the present experiments apparent that the formation 
cracks and the resulting disruption the lateral constraint the grains 
must make significant contribution the observed creep rate. The increasing 
creep rate with time the start the test would support the thesis that the 
tensile stress concentration across planes, such A-A Fig. tends 
increase with time. The formation crack would relieve the tensile stresses 
locally and thus change the stress distribution. This change stress distribu- 
tion would tend decrease the creep rate the grains involved. would 
appear that under the conditions the present experiments, during the first 
few minutes least, Andrade’s law for creep not obeyed. 

Steinemann (1954) gives evidence the recrystallization ice during the 
creep process. The ice used his experiments had granular rather than 
columnar structure. possible that for compressive stresses less than about 
kg/cm? temperature —10° creep and recrystallization the grains 
ensures that concentrations stress and imperfections will not develop the 
magnitude necessary for the initiation crack. 

Although the evidence indicates that the grain boundaries that control 
the formation the cracks and the direction the creep multigrained ice, 
the possible influence the easy slip direction, which tends parallel 
the direction the grain boundaries because the bias the crystallographic 
orientation, cannot ruled out. This bias may certainly have been con- 
tributing factor the lack creep the direction the long axis the 
grains. would also contribute the ability each grain deform under 
the constraints imposed its neighbors. 

Each time crack forms, the internal structure the ice altered. 
the present situation would appear that when enough cracks form, the 
ability the ice withstand shear stresses lowered the extent that 
large movements accompanied extensive crack formation occur along 
preferred planes which the shear stress maximum. When this condition 
reached the ice can considered have failed. The cracking associated 
with the planes shear failure responsible for the high cracking rate 
observed after the first hour for those test pieces that showed continually 
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increasing creep rate with time. probable that, the tests lower stress 
had been continued long enough, failure along the planes maximum shear 
would have developed these cases also, resulting sudden increase the 
cracking rate and corresponding increase the creep rate similar that 
observed concrete Riisch (1959). 

Because the movement that occurs the planes cracking Fig. 
and 7d, the creep measured extensometer similar that used would 
depend the location which mounted. obvious example this 
curve Fig. was also observed that the cracking tended occur 
one side the test piece, the creep would greater that side, that is, 
the test piece would develop bow. These factors are responsible part for 
the irregular dependence the creep the applied stress seen Fig. 
Even so, the observed dependence the creep rate the stress does agree 
well with the dependence found Steinemann (1954) for granular ice. 


CONCLUSIONS 


The purpose these observations was investigate, one temperature 
and under constant compressive load, the cracking that occurs ice during 
creep. particular, was desired investigate the sudden transition 
the curve that had been reported earlier and see this transition 
was some way related the internal cracking. From the present observa- 
tions the following conclusions are drawn: 

(1) During creep under uniaxial compressive stress greater than about 
kg/cm? applied perpendicular the grain boundaries ice —10° 
internal cracks form the ice. The rate formation these cracks depends 
the time and the magnitude the stress. For stresses below about 
kg/cm?, this rate goes through maximum within the first hour after the 
application the load and then decreases. The cracks tend long and 
narrow with their plane parallel the direction the grain boundaries and 
the applied stress. 

(2) The formation these cracks due tensile stress which developed 
perpendicular the applied stress because the relaxation shear stresses 
the grain boundaries. 

(3) The transient portion the curve due, part, the 
formation internal cracks. 

(4) The formation internal cracks weakens the structure until condi- 
tion attained where the specimen will fail along the planes maximum 
shear. 

(5) Because the logarithmic dependence the rate cracking 
the stress, this condition attained very rapidly stresses about 
kg/cm? and responsible for the sudden change the character the 
curve that observed this stress range. 

(6) Once failure begins along the planes maximum shear the cracking 
activity greatly increases, the cracks forming primarily the failure plane. 
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TABLES ROTATIONAL ENERGY LEVELS 
DEFORMED EVEN-EVEN 


ABSTRACT 


digital computer has been used extend calculations rotational energy 
levels deformed even-even nuclei treated the manner Davydov and 
Filippov. Tables are presented which may useful for comparison with experi- 
mental energy levels high spin such are realized Coulomb excitations 
heavy ions. 


Energy levels axially asymmetric nuclei, given 
Hamiltonian 


proposed Bohr (1952), have been examined Davydov and Filippov 
(1958), and Davydov and Rostovsky (1959). Here the mass transport 
parameter, and are shape parameters corresponding deformation 
order and the are operators the projections the total angular momen- 
tum along the body-fixed axes the nucleus. Certain energy levels are 
excluded symmetry conditions. 

previous paper, DeMille a/. (1959) calculated the energy levels for 
the then experimentally known excited states. important recent develop- 
ment Coulomb excitation techniques (Stephens, Diamond, and Perlman 
1959) has extended the observations and encouraged the present authors 
calculate the energy levels for states higher spin. The results such calcula- 
tions are presented Tables and II. These tables give the energy ratios 
E(n+)/E(2+) and also the energy the first excited state (in units 
that the ratios may converted energies desired. Figure 
shows the ratios function for the state and the 2’+ 
level. 

The calculations were performed with the aid IBM 650 digital com- 
puter. Calculated final results for the energy levels have accuracy 
least 0.02%. Spot checks were made for the levels higher spin using the 
tables King, Hainer, and Cross (1943), and the asymptotic method 
Golden (1948). 

Stephens, Diamond, and Perlman (1959), using multiple Coulomb excitation 
realized bombarding with ions, have succeeded establishing experi- 
respectively. Assuming these nuclei axially symmetric, they show that 
the level schemes may given 
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Energy ratios E(n+)/E(2+) and the energy E(2+) the "ground state rotational 
axially asymmetric even-even nuclei. given degrees and E(2+) units 
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TABLE I 


E(6+) E(8+) E(10+) 

oO 4.0000 3.3335 7.0000 12,000 18.333 26.000 

3 4.0220 §,3355 7.0000 12.000 18.333 26.000 

4 4.0392 3.3331 6.9987 11.996 18.323 25.979 

5 4.0613 3.3327 6.9964 11.989 18.306 25.944 

6 4.0887 3.3320 6.9929 11.978 18.280 25.891 

| 4.1212 3.3309 6.9869 11.959 18.235 25.799 

8 4.1590 3.3291 6.9774 11.929 18.163 25.654 

9 4.2022 3.3265 6.9633 11.885 18.059 25.446 

10 4.2510 3.32e7 6.9430 11.822 17.912 25.158 
11 4.3055 3.3175 6.9152 tt, 737 24.790 
12 4.3658 3.3104 6.8776 11.624 17.467 24.338 
13 4.4322 3.3010 6.8284 11.481 17.164 23.818 
14 4.5048 3.2887 6.7656 11.306 16.813 23.246 
15 4.5836 3.2729 6.6875 11,099 16.424 22.645 
16 4.6688 3.2529 6.5930 10.866 16.010 22.028 
17 4.7603 3.2280 6.4819 10.610 15.581 21.409 
18 4.8582 3.1975 6.3557 10,341 15,148 20.796 
19 4.9622 3.1608 6.2169 10.064 14,718 20.196 
20 5.0718 3.1277 6.0690 9.7869 14.296 19.613 
21 5.1864 3.0684 5.9164 9.5139 13.888 19.050 
22 5.3047 3.0137 5.7632 9.2499 13,498 18.513 
23 5.4251 2.9552 5.6138 8.9989 13.128 18.005 
24 5.5451 2.8950 5.4721 8.7650 12.785 17.535 
25 5.6613 2.8361 5.3420 8.5526 12.474 17.108 
26 5.7691 2.7814 5.2273 8.3666 12.202 16.735 
a7 5.8629 2.7342 5.1320 8.2127 11.977 16.426 
28 5.9365 2.6977 5.0601 8.0968 11.808 16.194 
29 5.9837 2.6746 5.0153 8.0246 11.703 16.049 
30 6.0000 2.6667 5.0000 8.0000 11.667 16.000 
E(14+) E(16+) E(18+) E(24+) 

35.000 45.333 57.000 70.000 84.333 100,000 

3 35.000 45.333 57.000 70.000 34.333 100.000 

4 34.962 45.270 56.899 69.847 84.111 99.687 

S$ 34.897 45.160 56.725 69.584 83.730 99.153 

6 34.300 44.996 56.467 69.198 83.175 98.382 

7 34.633 44.718 56.032 68.554 82,264 97.142 

8 34.373 44.290 55.376 67.605 80.958 95.420 

9 34.005 43.701 54.500 66.384 79.340 93.362 

10 33.513 42.940 53.415 64.928 77.481 91.075 
11 32.907 42.043 52.188 63.345 75.518 88.715 
12 32.200 41.042 50.868 61.686 73.502 86.322 
13 31.424 39.985 49.512 60.011 71.490 83.952 
14 30.606 38.903 48.147 58.344 69.499 81.612 
15 29.772 37.820 46.795 56.701 67.541 79.316 
16 28.937 36.747 45.463 55.087 65.620 77.063 
1g 28.110 35.691 44.156 53.504 63.737 74.855 
18 27.299 34.659 42.880 51.960 61.900 72.701 
19 26.508 33.655 41.637 50.457 60.112 70.604 
20 25.740 32.681 40.434 49.000 58.380 68.572 
2i 25.002 31.743 39.276 47.598 56.712 66.616 
22 24.296 30.849 38.170 46.260 55.119 64.747 
23 23.631 30.004 37.126 44.997 53.616 62.983 
24 23.013 29.221 36.158 43.824 52.220 61.344 
25 22.453 28.511 35.280 42.761 50.954 59.859 
26 21.964 27.890 34.512 41.831 49.847 58.560 
27 21.559 27.376 33.877 41.062 48.932 57.485 
21.254 26.990 40.484 48.243 56.676 
29 21.065 26.749 33.101 40.123 47.813 56.172 
30 21.000 26.667 33.000 40.000 47.667 56.000 
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TABLE 


Energy ratios the band axially asymmetric even-even nuclei. 


4 102.091 E( 2+) 104.425 407.697 106.091 409.364 108.092 
b | 65.162 67.496 259.982 69.162 261.649 71.166 
6 45.103 47.438 179.743 49.103 181.410 51.110 
10 15.937 18.285 63.079 19.937 64.749 21.998 
11 13.094 15.448 51.704 17.094 53.376 19.185 
12 10.933 13.297 43.060 14.933 44,734 17,065 
9.2543 11.631 36.340 38.017 15.438 
10.317 31.016 11.924 32.698 14.176 
15 6.8541 9.2683 26.729 10.854 28.416 13.190 
5.9811 8.4223 23.230 9.9811 24.924 12.420 
5.2608 7.7358 20.340 9.2608 22.043 
18 4.6608 7.1778 17.929 8.6608 19.643 11.359 
19 4.1572 6.7253 15.900 8.1572 17.629 11.010 
20 3.7320 6.3615 14.181 7.7320 15.928 10.751 
21 3.3717 6.0735 12.717 7.3717 14,487 10.565 
22 3.0659 : 5.8515 11.464 7.0659 13.263 10.437 
23 2.8068 5.6879 10.391 6.8068 12.227 10,353 
24 2.5888 5.5765 9.4725 6.5888 11.355 10.298 
25 2.4078 5.5128 8.6899 6.4078 10.631 10.257 
26 2.2610 5.4935 8.0305 6.2610 10.044 10.212 
27 2.1472 5.5141 7.4874 6.1472 9.5886 10,151 
28 2.0656 5.5677 7.0625 6.0656 9.2623 10,080 
29 2.0164 5.6335 6.7741 6.0164 9.0657 10,022 
30 2.0000 5.6667 6.6667 6.0000 9.0000 10.000 
3 730.517 1634.898 190.211 732.844 192.878 195.878 199.211 
4 411.364 916.819 110.424 413.697 113,095 116,090 119.433 
5 263.649 584.459 73.495 265.983 76.173 79.160 82.521 
6 183.411 403.922 53.435 185.745 56.124 59.099 62.485 
7 135.037 295.078 41.341 137.372 44.049 47.002 50.433 
8 103.648 224.449 33.492 105.983 36.230 39.150 42.652 
9 82.134 176.039 28.111 84.472 30.894 33.764 37.370 
10 66.753 141.426 24.263 69.094 27.110 29.908 33.659 
ll 55.383 115.832 21.416 57.726 24.349 27.050 30.995 
12 46.744 96.382 19.251 49.092 22.298 24.870 29.061 
13 40.032 81.261 17.566 42.386 20.758 23.164 27.655 
14 34.719 69.281 16.227 37.080 19,597 21.800 26.635 
15 30.445 59.635 15.146 32.816 18.724 20.685 25.899 
16 26.963 51.760 14.259 29.347 18.073 19.756 25.366 
17 24.095 45.254 13.521 26.494 17,593 18.966 24.974 
18 21.711 39.826 12.898 24.131 17.240 18,280 24.667 
19 19.718 35.257 12.367 22.162 16.980 17.671 24.398 
20 18.045 31.383 11.908 20.517 16.784 17.122 24.120 
21 16.642 28.081 11.507 19.144 16.623 16.618 23.794 
22 15.469 25.253 11.152 18.005 16.469 16.151 23.398 
23 14,502 22.827 10.836 17.069 16.295 15.717 22.942 
24 13.729 20.744 10.555 16.313 16.084 15.315 22.458 
25 13.150 18.960 10, 306 15.723 15.838 14.951 21.982 
26 12.771 17.444 10,090 15.281 15,582 14.632 21.549 
27 12.605 16.172 9.9129 14.976 15.347 14.367 21.182 
28 12.642 15.136 9.7790 14. 788 15.160 14,167 20.902 
29 12.837 14.355 9.6952 14.694 15.041 14.042 20.727 


30 13.000 14.000 9.6667 14.667 15,000 14,000 20.667 


| 
| 
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and give appropriate values for the constants and The 
nucleus has second known state character with energy 
790 kev (McGowan and Stelson 1957). The parameter may 
then given 10.0° the usual method establishing 
the ratio E(2’+)/E(2+), and the observed levels matched within few 
per cent. With this value corrections the above type are still necessary 
for closer agreement with experiment, but the constants and are about 
one-half the values required Stephens 

may significant that for and values may obtained 
from any two observed members the “ground state such that the 
remaining levels this band may found within the experimental accuracy, 
without assuming any corrections. DeMille a/. (1959) have shown that this 
true for many deformed nuclei. The values obtained are commonly 
few degrees larger than those given the ratio 
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CENTRAL THREE-BODY FORCES HEAVY 


SMITH? AND ROBERT SHARP 


ABSTRACT 


Three-nucleon potentials are calculated using the renormalized, static theory 
Chew and Low. These potentials are used evaluate the three- body contribu- 
tion the total energy per nucleon nuclear matter function nuclear 
density. found that while the three-body energy greater than that predicted 
unrenormalized theories about one order magnitude result 
multiple-scattering effects, its dependence the nuclear density the region 
the equilibrium density very weak. Three-body forces are therefore not 
expected change the saturation properties nuclear matter predicted 
hard core potential any appreciable extent. 


INTRODUCTION 


Three-body forces have been interest nuclear physics for considerable 
time (Janossy 1939; Primakoff and Holstein 1939; Drell and Huang 1953; 
Brueckner, Levinson, and Mahmoud 1954; Fujita and Miyazawa 1957). 
first was believed that many-body forces might needed account for 
the saturation nuclear densities. More recently has been shown, however, 
that the inclusion phenomenological infinite repulsive core the two- 
body potential sufficient produce equilibrium nuclear density which 
reasonable agreement with experiment (see, for example, Gomes, Walecka, 
and Weisskopf (1958)). However, still interest know three-body 
forces are important, and so, whether agreement with experiment improved 
taking them into account. 

Two-body nuclear forces calculated with unrenormalized, pseudoscalar 
meson theory and fixed, extended sources were studied extensively 
Brueckner and Watson (1953). The three-body force obtained from this theory 
was calculated Brueckner, Levinson, and Mahmoud (1954). the non- 
relativistic approximation, where the effects pair 
formation are small, the three-body forces are found negligible; fact, 
using the approximation methods Drell and Huang (1953) (which involve 
the neglect Coulomb, nuclear surface, and core correlation effects) 
evaluate the expectation value the potential, they are found produce 
only about 0.1 Mev repulsion per particle. 

That these estimates the three-body energy are not reliable was pointed 
out Chew (1954) and Klein (1957). The reason that processes which 
involve the exchange two more pions are enhanced multiple scattering 
resonance will occur when the energy one the virtual mesons present 
the intermediate state equal the energy the exchanged meson. This 
the so-called (3,3) resonance—the resonance which occurs the 
nucleon channel with total spin 3/2 and total isotopic spin 3/2. 
Klein pointed out, the actual scattering this channel greater than the 
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Born approximation factor more. therefore seems likely that 
the three-body energy per nucleon nuclear matter, with multiple scattering 
taken into account, will appreciably greater than that predicted 
unrenormalized theory (i.e., when multiple scattering neglected). Thus, 
calculation the three-body energy per nucleon using the renormalized, static 
theory Chew (1954), and Chew and Low (1956), which takes these multiple- 
scattering effects into account, appears justified. 


CALCULATION THE POTENTIAL 
order calculate the three-body potential lowest order (two pions 
exchanged), the following type Hamiltonian chosen: 


(2) 


This natural extension the Hamiltonian used Sharp (1958) for the 
two-nucleon potential the problem three nucleons. The first six terms 
(1), which are taken the unperturbed Hamiltonian, comprise the 
Hamiltonians for three independent pion fields coupled nucleons and 
respectively; for example, H}, the energy the pions associated with 
nucleon and couples these pions nucleon The final six terms (1) 
give the interactions between the three nucleons; that is, allows nucleon 
absorb pion which has been emitted nucleon and on. These terms 
are treated perturbation. this means, the ‘‘renormalization 
the nucleons are treated exactly because they are contained the unperturbed 
Hamiltonian, and only the interaction between nucleons treated 
perturbation. The nucleon energies not appear because, the static 
approximation, the zero energy can chosen the nucleon rest energy. 
The interaction Hamiltonians are defined the following expressions 


(4) 


where the momentum the pion and charge state, normaliza- 

the energy the pion, and the momentum cutoff 
factor which gives rise extended nucleon sources for the pions. the 
unrenormalized coupling constant, and the annihilation operator for 
pion type the spin and isotopic spin operators respectively 
for nucleon 

The unperturbed, Hermitean part the Hamiltonian gives complete set 
and on. The unperturbed ground state this represents three 
nucleons, and free pions. (As used here, pion one which 
exchanged between nucleons.) 
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The potential due the exchange two mesons now given directly 
conventional second-order perturbation theory: 


where defined (2). The exchange one meson does not contribute 
the three-body force, and processes involving the exchange more than two 
mesons will neglected this calculation. The prime the sum indicates 
that least one m2, must state with one more pions present; 
that is, the unperturbed must not included. 

order simplify the calculations, choose nucleon that with two 
vertices, and sum over all permutations the three nucleons after the 
potential has been calculated. 

now invoke the one-meson approximation Chew and Low (1956). 
Thus, retain our potential only those intermediate states m2, m3) 
which one and one-meson state, and the others are real nucleon 
states. There are four essentially different Feynman diagrams which contribute 
our potential, plus number others obtained permuting the three 
nucleons; the basic diagrams are shown Fig. The term the perturba- 


3 


Fic. Diagrams the various processes contributing the three-body potential. The 
terms the Hamiltonian causing the absorption the mesons are written above. All other 
diagrams obtained permuting the three nucleons must also included. 


tion Hamiltonian which causes the absorption the pion written above 
the pertinent pion line (its emission accomplished the unperturbed 
The three-body potential can now written down 
inspection; 


Ws 


plus all terms obtained permuting the three nucleons. |s) one-meson 
state nucleon with momentum and the perturbation Hamiltonian has 
been written using equation (3). 
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order evaluate this expression for the potential, use the following 
relationships for the matrix elements. They are rearrangements equations 
(2.15), (5.32), (5.35), (5.37), (6.6), and (6.13), the article Wick (1955). 


where the renormalized coupling constant, the solution the Low 
(6.12). Sums are understood over all repeated indices, and over all When the 
sums over momenta are replaced integrals, 


the integrations over angles can performed follows: 


where 
and sin 


These results are applied the equation (6), resulting the 
following expression: 


q if 
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and 


the factors parentheses, the Latin indices refer spin, the Greek 
isotopic spin; before, sums are understood over all repeated indices. The 
unit vectors are now relative co-ordinates defined 


now make use the effective range approximation Chew, Goldberger, 
Low, and Nambu (1957). This amounts the assumption that the (3,3) 
resonance dominates the scattering picture, that the sum over the 
above equations, only contributes significantly. The potential can 
simplified the use the integral equation Low, which is, for 


(8) 


where defined Wick (1955). The second, term under 
the integral small because and contain resonances, and 
small; also, the energy denominator large, that this term can 
neglected. the limit large the Low equation becomes 

2 2 2 

where the resonance energy (see Chew, Goldberger, Low, and 
Nambu 1957). 

order simplify the potential, the Low equation the form equation 
(8) (with the crossing term omitted) substituted into the first two terms 
equation (7). found that the resonant energy denominators the fifth 
and sixth terms (7) are just cancelled similar parts the first two terms. 
Upon application (8a) the resulting expressions, the potential found 
assume the form 


where has been neglected comparison 


(8a) 
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Tensor forces are not expected produce any appreciable effect heavy 
nuclei, because adjacent low-lying states which might have been admixed 
are already occupied. fact, the purely two-body problem which concerned 
them, Gomes, Walecka, and Weisskopf (1958) found that the neglect tensor 
forces underestimates the binding energy roughly 15% the 
potential energy; for our purposes, this quite negligible. therefore average 
over the angular dependence the potential, which contained the factors 


(9) 


plus imaginary terms which are completely antisymmetric exchange 
particles and and vanish when all possible nucleon permutations are 
taken into account. 

The potential therefore 


Mev. Gaussian form cutoff was chosen, and these integrals were 
performed numerically. 


THREE-NUCLEON CHARGE-SPIN STATES 


There are various prescriptions for writing down states describing the 
different ways which the spins and isotopic spins three interacting nucleons 
‘can combine (see, for example, Verde (1957) and Sachs (1953)). propose 
new representation, using spin states which are eigenfunctions the total spin 
the component spin Sz, and the operator and isotopic spin 
states which are eigenfunctions the analogous isotopic spin operators. This 
representation has the advantage simplicity when used with the Hamiltonian 
equation (1), well the fact that the states are eigenfunctions the 
cyclic three-body exchange operator. 
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The spin states are 


123 
X33 


where and are cube roots unity. The subscripts are and respectively 
the superscript indicates the fact that the eigenvalue when 
operating the 3/2 state vectors zero, while indicates that its 
eigenvalues when operating the 1/2 state vectors are The 
isotopic spin states, which are derived analogous manner, will labelled 
spin states, obtained combining the charge states and spin states all 
possible ways. choose the following combinations: 


(1) V2 X1, xi, 1, (4) 
(3) (16) 
(5) xt att Tz (8) 


The states labelled (1) are completely antisymmetric upon interchange any 
two particles, while those labelled (2) and (3) are completely symmetric. 
Those labelled (4) and (5) are states mixed symmetry. The number such 
states shown parentheses the right; these numbers are measures the 
statistical weights the corresponding states. 
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The expectation value the potential operator these various states must 
now evaluated. Since all the spin dependence the 
potential contained the factor equation (9), will sufficient 
simply take matrix elements this operator. These are shown Table 


TABLE 
Matrix elements operators 


The entry the column labelled reality 


Thus, can diagonalized simultaneously with the Hamiltonian. 

The three-body force, because its short range, should strongest 
symmetric spatial states, that is, antisymmetric states. But 
seen from Table that the expectation value the antisymmetric 
states vanishes, that what might expected the dominant 
three-body force fact zero. must therefore consider contributions from 
states which are not completely antisymmetric upon exchange 
any two nucleons. 


IV. THE SPATIAL WAVE FUNCTION 


Besides the spatial states which are symmetric interchange any pair 
nucleons, which have been found not contribute the three-body force 
due the vanishing the charge-spin matrix element associated with them, 
there exist spatial wave functions which are antisymmetric interchange 
any pair particles, and wave functions mixed symmetry. The antisym- 
metric spatial states will not contribute significantly the three-body force 
because, roughly speaking, the antisymmetry the wave function keeps the 
nucleons outside the range the three-body force. This can demonstrated 
constructing function which expresses the correlations imposed the 
exclusion principle; the correlation function derived manner which 
analogous that Wigner and Seitz (1933) for two-body correlations due 
the exclusion principle. 

The spatial wave function chosen 


mn 


where the Levi-Civita tensor density. The one-particle wave functions 
are chosen plane waves because real scattering can take place 
nuclear matter; all adjacent states are occupied. averaging over all 
momenta from zero the Fermi momentum one obtains the correlation 


function 
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i= 
where 
3(sin cos 


This function shown graphically Fig. function for various 
values assuming two configurations: (i) the particles are equidistant 


Fic. The function which expresses correlations due the exclusion principle the anti- 
symmetric spatial states plotted against Curve (a): three colinear nucleons, 1.69 
curves (c), and (d): three equidistant nucleons, 1.69, 1.48, and 1.27 
respectively. 


from each other (7; and (ii) the particles 
(r1 = 2re = 2r3). 

the expression for the potential (equation (10)) that the term 
with the longest range drops half its maximum value 
cm, which point the correlation function P(r) for 
(the case when P(r) largest) only 0.12, and rapidly becomes smaller 
closer distances. There therefore very little likelihood the nucleons 
getting close enough together within the range the three-body force. 
Thus, consistent with previous approximations neglect the effect 
antisymmetric spatial states. 

seen therefore that the dominant three-body force arises the spatial 
states mixed symmetry. order calculate the binding energy due 
this force, function which takes into account the correlations induced the 
wave function the exclusion principle must again calculated. Again using 
the methods Wigner and Seitz (1933), the correlation function found 


P(r) 


where has been defined previously. seen that this correlation function 
goes zero only when all three the r’s vanish. Thus, the wave function 
vanishes only linearly small distances instead cubically, was the case 
for completely antisymmetric spatial wave function. 


3 3 
1.0 
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Account must also taken the correlation induced the spatial wave 
function the infinite repulsive core. High energy scattering data indicate 
infinite inside the core, the wave function must vanish this region. 
real scattering allowed nuclear matter, however, due the effect the 
exclusion principle. Thus, there can permanent phase shift large 
distances; the real wave function must approach the free-particle wave 
function asymptotically. The distance which the difference between the 
real wave function and the free-particle wave function becomes small will 
called the accordance with the nomenclature 
Gomes, Walecka, and Weisskopf (1958). This distance will order 

order introduce this correlation, will sufficient for our purpose 
multiply the previously calculated function which approaches 
unity asymptotically, but which close unity distance order 


TH—Te 


This function plotted Fig. for various values 


1.0 


(r) 


r 
Fic. The function which expresses correlations due the infinite repulsive core plotted 


The square the spatial wave function now 


CALCULATION THE THREE-BODY BINDING ENERGY 


order calculate the three-body contribution the binding energy, 
the expectation value the potential evaluated the previously specified 
three-body states nuclear matter. Thus, the binding energy per particle 


2 
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The which used here the potential energy before summing over 
all permutations the three nucleons, given equation (10). The sum over 
all nucleon permutations introduced the factor 
one factor then divided off obtain the binding energy per particle. 
The factor where the nuclear volume, used normalize the 
one-particle wave functions nuclear volume rather than unit volume. 
Table summarizes the results this calculation. 


TABLE 
Nuclear density relationships 
E/A 
Mev per particle 

1.27 0.138 1.20 1.24 
1.48 0.219 1.03 1.49 
1.69 0.328 0.90 2.82 


Note: ro = RA~'3, where R is the nuclear radius. 


The binding energy plotted function Fig. along with the 
total two-body energy per nucleon calculated Gomes, Walecka, and 
Weisskopf (1958) the Born approximation their 
and the sum these two energies. Due the weak dependence 
the three-body energy the region the minimum the two-body 
total energy curve, the equilibrium nuclear density not changed any 
appreciable extent. can concluded from this that the saturation nuclear 
density due principally the hard core, since the three-body force (which 
the strongest many-body force) has very little effect the equilibrium 


density. 


VI. CONCLUSIONS 


The aim the present work has been twofold: firstly, construct three- 
body potential using coupling static nucleons the meson 
field; and secondly, using this potential, calculate the effect three-body 
forces the saturation properties nuclear matter, assuming two-body 
potential with infinite repulsive core. has been found that the three- 
body force has moderately short range, would expected from the fact 
that the process involves the exchange two pions. 

new method specifying the states three nucleons has 
been introduced. This representation has the advantage over previous schemes 
that the symmetry characteristics the states are immediately apparent, 
well the fact that they are eigenstates the Hamiltonian used this 
calculation. Antisymmetric states comprise 1/16th, symmetric 
states 5/16ths, and states mixed symmetry the remaining the 
total states. 

was shown that antisymmetric nuclear three-body spatial states could 
neglected, because the antisymmetry the spatial wave function keeps 
the nucleons outside the range the three-body force. The contribution 


| 
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(Mev) 


Fic. The binding energy per nucleon function the Fermi momentum Curve (a): 
three-body binding energy; curve two-body binding energy calculated Gomes, Walecka, 
and Weisskopf (1958); curve (c): the sum these energies. 


the symmetric spatial states, which might expected account for the 
dominant part the binding, was found exactly zero, due the fortui- 
tous vanishing the spin isotopic spin matrix elements the potential 
these states. Thus, only the spatial states mixed symmetry contribute 
significantly the binding energy. these states, the three-body force 
repulsive, agreement with unrenormalized calculations. 

The belief that multiple-scattering effects, which give rise the resonance 
the (3,3) scattering channel, increase the magnitude the three-body force 
borne out this calculation. The three-body repulsion obtained here 
greater about one order magnitude than that obtained Brueckner, 
Levinson, and Mahmoud (1954), who neglected all ‘‘renormalization 
including multiple scattering. 

Brueckner and Watson have suggested that iterated second-order terms, 
which appear implicitly the potential obtained here, should not included 
the fourth order calculation. While agreement this subject has not 
yet been reached, would simple matter calculate the contributions 
these processes the potential standard perturbation methods, and 
subtract them from the potential which has been calculated. 

The three-body potential resulting from the exchange three mesons 
‘cannot calculated present-day techniques. The difficulties arise 
applying the one-meson approximation Chew and Low those terms 
whose diagrams have three vertices one nucleon. One possible means 
including these terms would the use experimental cross sections for 
multiple meson production place more fundamental calculation the 
nucleon three-meson interaction. 

While the neglect angle-dependent forces seems reasonable, investiga- 
tion into their exact effect might prove fruitful, because the dependence 


3 
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these forces the nuclear density may quite different from that the 
central forces. Since nuclear recoil has also been neglected, velocity- 
dependent forces arise. means certain that these are not required; 
they may very well important both two- and three-body interactions. 

can concluded that central three-body forces play relatively small 
role the saturation nuclear densities, because the three-body binding 
energy per nucleon has little effect the equilibrium nuclear density obtained 
from two-body forces. The repulsive core remains the dominant factor 
explaining the saturation properties nuclear binding energies. 
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APPENDIX 


The expression for the three-body binding energy (equation (11)) 
integral over the spaces the vectors and r3. complicated the fact 
that appear correlation functions depending the vector which 


(12) 


where the cosine the angle between and The integrations over the 
two spaces are therefore not separable. was accordingly decided replace 
the and equation (12) effective values, order make the integrals 
separable. different effective was calculated for each the different 
products appearing the expression for the binding energy, for each value 
ky, calculating the expectation value for the various terms the 
potential. 
When this done, equation (11), which the form 
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becomes 


where the functional dependence and has been replaced the 
effective values determined above, thereby making the integral separable. 

The integrations over and extend from the core radius infinity. 
However, the integration over from does not prevent from 
becoming smaller than the core radius That this does not alter the value 
the integration appreciably can seen follows. The core will make itself 
felt during relatively small proportion the time, because can become 
small only when and are approximately equal magnitude. The worst 
case arises when both and are equal the smallest effective distance 
calculated above; this case, found that only 7.5% the range 
integration should omitted and are exactly equal magnitude. 
This therefore represents upper limit the error. 

The integrals themselves were performed 650 digital computer. 


THE TEMPERATURE DEPENDENCE 
THE SHAPE PARAMAGNETIC RESONANCE LINES! 


McMILLAN? AND OPECHOWSKI 


ABSTRACT 


theory the temperature dependence the shape paramagnetic resonance 
lines developed for temperatures sufficiently low make the effects the 
spin-lattice interaction and thermal expansion negligible. General expressions 
for the first and second moments the line shape function functions the 
temperature are obtained, and the approximations used are discussed detail. 
These expressions are applied the case effective spin and 
illustration the general theory numerical results are given for the paramagnetic 
resonance lines nickel fluosilicate. this case the moments become strongly 
temperature dependent below 10° 


INTRODUCTION 


sufficiently low temperatures the effect the variation the tempera- 
ture the paramagnetic resonance absorption mainly due the tempera- 
ture dependence the average occupation numbers (as given the Boltzmann 
factors) the various energy levels the absorbing system. The effects 
connected with lattice vibrations, the and thermal 
expansion, become small these temperatures. this paper neglect them 
altogether and treat the lattice perfectly rigid. other words regard 
the position and shape paramagnetic resonance line determined solely 
the nature the paramagnetic unit systems (‘‘spins’’) which constitute 
the paramagnetic crystal, and the weak interactions these unit systems 
this assumption, derive expressions which 
give the temperature dependence the first and second moments the 
absorption line for the case which the paramagnetic ion subjected 
arbitrary crystalline electrostatic field, and external constant 
magnetic field. 

The temperature dependence the first two moments line has been 
considered Pryce and Stevens (1950), and Kambe and Usui (1952) 
Pryce and Stevens have derived some general approximate expressions con- 
taining Boltzmann factors but when specializing these expressions they have 
considered the temperature infinite. Kambe and Usui have set quite 
general exact expressions for the moments, which shall briefly rederive 
Section and use the starting point our calculations. However, Kambe 
and Usui have applied their general formulae only the case which 
crystalline electrostatic field present. From their results is, however, quite 
clear that the well-known Van Vleck (1948) temperature-independent formulae 
for the moments could longer used interpret experimental data 
paramagnetic resonance lines temperatures the order 10° lower, 
should such data sufficient precision become available. 
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The results the present paper can thus regarded straightforward 
extension the results Pryce and Stevens, and Kambe and Usui. fact, 
use the methods calculation developed these authors, and not 
claim have introduced any essentially new However, discuss 
more explicitly than previous authors the meaning various approximations, 
and especially the temperature-independent approximation for the 
moments. 

the calculations which follow assume that the amplitude microwave 
radiation very small. other words, assume the absence saturation 
effects. This may seem incompatible with the assumption that the 
lattice interaction negligible which implies long relaxation 
time and, consequently, accumulation energy the spin system. How- 
ever, the relaxation time may very well short enough prevent 
observable saturation effects occur and still long enough not contribute 
the shape the absorption line appreciable way. 

The case nickel fluosilicate treated Sections and this paper 
intended mainly illustration the application the very complicated 
general formulae given Section 


EXPRESSIONS FOR THE FIRST AND SECOND MOMENTS THE 
LINE SHAPE FUNCTION 

implied what have said the Introduction, use this 
paper the customary method moments (Waller 1932; Broer 1943; Van Vleck 
1948) order characterize the position and shape paramagnetic reso- 
nance line. 

shall define this section shape such that f(v)is 
proportional the imaginary part the high frequency susceptibility 
paper expressions the type 


order able define the line shape function must first describe 
more explicitly the system” for which our formulae will valid. The 
lattice, each spin being defined some equiva- 
lent expression. The spin-Hamiltonian incorporates, usual, the effect 
the external constant magnetic field and the crystalline electrostatic field. 


The Hamiltonian the spin system then given 


where 


N 


9g” = 9° 


i=l 


and describes the interactions (which include magnetic dipole 


> 
> 
> 
> 
» 
> 
’ 
i 
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and exchange interactions), and describes the effect the high frequency 
oscillating magnetic field the spins. shall assume that the interactions 
are weak, i.e., 


and also 


that the Hamiltonian the spin system can written approximately 


shall consider the unperturbed energy and the perturbation 
energy. shall treat causing transitions between energy levels 
with two qualifications however. Firstly, according the well- 
known argument Van Vleck (1948), which will not repeat here, the 
term the Hamiltonian must appropriately before 
can use for calculating the moments the line shape function f(v). The 
truncation necessary order eliminate the weak ‘‘secondary”’ 
lines from and consists discarding that part which does not 
commute with shall write for the truncated Secondly, 
further truncation necessary, pointed out Kambe and Usui (1952), 
when field present, namely truncation This truncation 
has with the fact that when crystalline field present the energy levels 
unperturbed spin are, general, not equidistant, and, hence, several 
absorption lines comparable intensity, that is, several lines, will 
observed. order calculate the characteristics for only one these 
lines, one has then modify way make cause transitions 
between only two levels spin. This kind truncation will denoted 
circumflex over (or over the operators which enter the expression for 
shall write shall define more explicitly later on, 
terms its matrix elements. 

Before proceeding, must say few words about the notation used 
for the Boltzmann factors. Since over 100 Boltzmann factors occur the 
equations the present paper, will consistently use the abbreviations 


—1/kT 


expression, may convenient use the notation for analogy 
with the notation for the exponential function, i.e. define function 


bol follows: 


shall label the eigenvalues the Hamiltonian the spin system 
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where the amplitude the oscillating magnetic field and the 
component the magnetic moment the spin system the direction the 
oscillating field. Then, thermal equilibrium, the power absorbed 
the spin system from the oscillating field will proportional 


nn’ 
where means summation over all states and |n’) such that 
nn! 


the other hand, have the well-known relation between and the 
imaginary part the high frequency susceptibility: 


P(v) 


now define the line shape function f(v) times normalization 
constant, 


0 


then this equation combined with the two immediately preceding equations 
gives 


where means summation over all states such that 


n’ 


Hence, using definition (2.1), can write the expressions for the first two 
moments follows: 


shall calculate these two moments. usual this kind calculation, 
shall first rewrite (2.2) and (2.3) invariant form, involving traces 
operators. Following Kambe and Usui (1952), get rid the rather 
awkward summation introducing the operators and which 
are defined 


y 
ly 
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otherwise, 
otherwise. 


Equations (2.2) and (2.3) can then rewritten, has been shown Kambe 
and Usui, the following invariant form: 


(2.4) 
trace 


trace 


Equation (2.4) can also written 


trace 


trace 


should noted that equations (2.4) and (2.5) have been derived without 
resorting any assumptions the temperature These equations are 
quite general and are valid all temperatures. Furthermore, they could 
applied not only the case electron paramagnetic resonance but also 
the case nuclear magnetic resonance. However, the latter case the region 
temperatures for which the usual temperature-independent formulae 
Van Vleck cease valid far below and hence probably 
immediate experimental interest. 

shall now rewrite equations (2.4) and (2.5) terms and rather 
than and M_. this end shall introduce certain projection 
operators the way Pryce and Stevens (1950) have done before. 

define projection operator, such that 


where satisfies the eigenvalue equation 


consequence (2.6) have that 
(2:2) Py Py Py Oya. 


Now clear that any operator can written 


2 Q 


w=1 n=l 
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where the number different eigenvalues Using this and the fact 


(2.8) 


Let now consider shall suppose that equations (2.4) and (2.5) 
give the first and second moments the primary line the frequency 
say. That is, consider the energy difference the two levels 
the unperturbed spin between which induces transitions. The circumflex 
has been introduced over order eliminate from the expressions for the 
moments any the primary lines which are not interested. Thus 


have 
Then can write 


where means summation over all values and wherein Eg— 
Let now suppose that also satisfies the equation 


(This possible since 0.) Then, not too large, have 
pairs states this the case then can write 


a.B 


Using equations (2.7) (2.11), equations (2.4) and (2.5) can rewritten, 
after some straightforward rearrangements, follows 


where 


trace 


a,p 


; 
| 
| 
| 
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now define 
That is, denotes the /th moment f(v) about the frequency should 
noticed that not general equal the /th central moment 


this being defined 
Using equation (2.14) and the fact that 


equations (2.12) and (2.13) give 


where the expressions for and differ from those for and only 
the replacement everywhere except the exponentials. 


DISCUSSION VARIOUS APPROXIMATIONS INVOLVING 
THE TEMPERATURE 

The expressions for the moments discussed the previous section not 
involve any assumption about the magnitude the temperature. 

well known, most physical situations experimental interest the 
temperature such expressions can regarded infinite. order obtain 
the limiting values these expressions for not sufficient replace 
the exponentials unity, because this would lead only indeterminate expres- 
sions the form 0/0 glance equations (2.4) and (2.5) shows that this 
so, since the trace commutator two finite matrices always vanishes). 
However, expanding the exponentials into powers one easily finds 
that the correct limiting values arise from terms linear and are 


given 


where 


(3.5) 


and 


P 
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a8 


not difficult show that these last three equations can obtained 
from 


(En —E,)'|(n| |n’)|? 


then using equations (2.7) (2.11) and equation (2.14). the Appendix 
show how can arrive equation (3.6) from equation (2.1). 

Equation (3.6) obvious generalization Van Vleck’s (1948) equation 
(3) and for this reason call the Vleck” approximation for the 
moment. 

for odd, often the case (see Van Vleck (1948)), then, 
equations (3.1) and (3.2) can written 


lim (Av) 


lim 
exact limiting value (Av) not zero predicted the Van Vleck 


Tc T3020 
usually done. 

The approximation consisting treating the temperature infinite can 
expected reasonably good for most paramagnetic salts down 
temperature few degrees Kelvin. breaks down soon 
becomes comparable with unity. few additional remarks about this point 
will made the end Section 4.) Hence, much better approximation 
will consist making use the fact that the interactions are weak, i.e., 
and expanding the exponentials equations (2.16) and (2.17) 


follows: 


(1) 


obtaining (3.7) use has been made equation (2.15). shall call the 
approximation consisting keeping only the first term the square bracket 
“the first approximation”, and the approximation consisting keeping only 
the first and the second terms the bracket second 
straightforward calculation shows that the exact equations (2.16) and (2.17) 
reduce, the first approximation, 


¥ 
n n’ 
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where 


trace 


Equations (3.8) and (3.9) are exactly equivalent equations (1), (2), and (3) 
Pryce and Stevens’ paper (1950). 

may pointed out that letting infinity equations (3.8) 
and (3.9) one obtains the Van Vleck approximation for the moments, i.e., 
equations (3.3) and (3.4) rather than equations (3.1) and (3.2). The latter 
equations contain small additional terms arising from the here neglected 
second term the bracket equation (3.7). 


EXPRESSIONS FOR THE MOMENTS THE FIRST APPROXIMATION 
FOR SYSTEM IDENTICAL “SPINS” 

shall now apply equations (3.8) and (3.9) the case when the unper- 
turbed spin has energy values which are all non-degenerate. 
shall suppose that pairs these energy values have energy separation 
hv*. Now, let suppose the pairs levels have energies 

shall denote the set positive integers (x, and have 

order evaluate the traces and perform the summations equations 
(3.8) and (3.9) use the same methods used Pryce and Stevens (1950) 
Section their paper. That is, take 


N 


2 
then evaluate the traces with respect the functions 
where, for example, satisfies the eigenvalue equation 


Any eigenvalue can written 


R 
>on, a, 


r=1 


where 


R 
and 


r=1 
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that the summations over the and can replaced summations 
over 
the method not new, and the calculations are rather lengthy and 
somewhat cumbersome, give here only the results.* 
introduce the following notation: 
denotes the single spin state having the energy 
denotes the matrix element the single spin operator between 
the single spin states |e) and |e’); 
denotes the matrix element between the state 
the spins and and the state the same pair spins; 


means summation over all states |e) single spin, 


means summation over all integers the set defined the beginning 


this section, 


means summation over all positive integers and such that 
ab 


means summation over all positive integers and 
ab 


means summation over all positive integers and such that 
ab 


means summation over all and independently, excluding those terms 


means summation over all and independently, excluding those terms 
ijk 

for which any two the three summation indices are equal. 

Further, 


= 


With this notation the first two moments are: 


e f 

e 


*For some details the calculations one may consult McMillan (1959). 


q 


| 
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where 


Ais(ef) 


have verified that the absence crystalline field, i.e. the case which 
each unperturbed spin has equidistant energy levels, equations (4.1) 
and (4.2) become identical with the expressions given for this case Kambe 
and Usui (1952). 

pointed out Section equations for (Av), and reduce, for 
the formulae for the moments the Van Vleck approximation (equations 
(3.3) and (3.4)). glance equations (4.1) and (4.2) makes clear that 
can regarded soon the energy differences are 
negligible when compared k7. more careful analysis the expression 
(4.2) shows that 


the smallest and the largest single spin energy value. 


APPLICATION THE FIRST APPROXIMATION FORMULAE FOR THE 
MOMENTS 

General equations (4.1) and (4.2) for the first and second moments are 
complicated, and not easy see from them glance what sort 
temperature dependence (Av); and they imply. 

will specialize equations (4.1) and (4.2) the case which each un- 
perturbed spin has not equidistant energy levels, and the set defined 
Section consists one element. For convenience, will label the energy 

Equation (4.1) gives then the following formula the case absorption 
line from the ground state the first excited state the unperturbed spin: 
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Similarly one obtains from equation (4.2) expression for the second 
moment the same line. not give this expression here because not 
much simpler than equation (4.2) itself. write later for two special 
cases (see equations (5.9) and (5.10)). 

further specialize the problem assuming that the spin-Hamiltonian 
single spin has axial symmetry, that the constant magnetic field 
directed along the axis symmetry which choose the and that 
the eigenstates the spin-Hamiltonian are also eigenstates the operator 
corresponding the z-component the spin. assume that the 
interaction sum magnetic dipole interaction and isotropic exchange 
interaction. can then write the interaction operator 
its usual form: 


the direction cosine the line length joining spin spin relative 
the z-axis; are the exchange interaction constants, which shall assume 
all equal spins and are nearest neighbors, and zero otherwise; 
and are the g-factors perpendicular and parallel the symmetry axis; 


the Bohr magneton. 


Substituting the appropriate matrix elements (5.4) equation (5.1) 


one obtains the case the effective (R=2 


(the terms involving cancel out), and the case the effective spin 
S=1 


For the second moment the same line one obtains the case 


| 
3 
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tj 


tj 
+2(b% (b%—b"*)A 17D jx]. 


the shape the sample spherical, and the lattice spins cubic 
(each spin having nearest neighbors) one has 


Equations (5.7) and (5.8) become then 


obtain some still more specific expressions for the moments consider, 
example, the case nickel fluosilicate. this paramagnetic crystal there 
only one Ni** ion the unit cell, which greatly simplifies the problem. The 
paramagnetic resonance absorption lines nickel fluosilicate have been 
observed Holden, Kittel, and Yager (1949), and Penrose and Stevens 
(1950), and the theoretical expressions for the moments have been obtained 
Ishiguro, Kambe, and Usui (1951), using the Van Vleck approximation. 

While Ishiguro, Kambe, and Usui assumed exchange interaction nickel 
fluosilicate isotropic, would appear from Ollom and Van 
(1951) theoretical analysis experimental evidence that fact strongly 
anisotropic. some points that analysis were rather uncertain, have 
assumed for simplicity that exchange interaction isotropic, and have 
used equations (5.3) (5.6) they stand. However, after our calculations 
were essentially concluded new experimental evidence support the 
anisotropy exchange interaction nickel fluosilicate was obtained Walsh 
(1959). seems thus that our quantitative predictions concerning nickel 
fluosilicate ought not taken too literally; but they should help any 
case clarify the situation this point when relevant experimental data 
become available. 

has long been known (Becquerel and Opechowski 1939; Penrose and 
Stevens 1950) that each ion nickel fluosilicate has effective spin 
and that the spin-Hamiltonian ions can written 


(5.14) s” = { —Blg. S,iHy) +g) 


where and are positive constants. consider here only the case 
which the constant magnetic field parallel the symmetry axis the 
crystal that 

(0) 

(5.15) 
and the three eigenvalues are 
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where, since have introduced the convention a3, assumed 
that large enough make There are only two absorption 
lines, corresponding 


and 


shall assume, has been done Ishiguro, Kambe, and Usui (1951), 
that the lattice the ions can approximated simple cubic lattice 
(lattice constant d), and shall take for granted the following numerical 
result due them: 


ij 


shall further assume that the crystal has the shape sphere that 
equations (5.11) hold. From these equations one easily derives the following 
identities: 


tye 
these assumptions the first two moments the line corresponding 
obtained from equations (5.8) and (5.10), are 


The first two moments the line corresponding can 
found appropriately relabelling the eigenstates equation (5.1) and 
the equation (not written here explicitly) for the second moment: one has 
replace and the matrix elements and Boltzmann factors. 
One then obtains using the same assumptions before 


(5.21) 
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The first and second moments for the two lines, given equations 
(5.19) (5.22), are plotted function the temperature Figs. 


300 Approx 
Approx 


gauss 


ISt Approx 


200 
| 
| 


the temperature the first and second approximation for the positive and negative sign 


the exchange interaction constant Conversion factor: gauss corresponds 2.19 
Fic. The first moment the line corresponding asa function 
the temperature. For further explanations, see the caption Fig. 


The numerical values these constants are taken from Penrose and Stevens 
(1950). fact temperature dependent between room and liquid hydrogen 
temperatures, but below 20° seems become constant, and this 
constant value that have used. regard and also temperature 
independent, which probably they really are the very low temperatures 
are interested in. Their numerical values are taken from Ishiguro, Kambe, 
and Usui (1951). They obtained two different values for the two 
lines using the experimental data Holden, Kittel, and Yager (1949): 
|4| erg for the line corresponding have 
used these values our numerical calculations. From the point view 
our theory would have been more logical course use the same value 
for both lines. 

Figs. and have also shown the second central moment 

The conversion factor from gauss erg has been chosen equal i.e. 

Some rather obvious comments Figs. are made Section 


Temperature 
300 
400« . 
— 2 
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70; 


Temperature 


600 

500 

400 

300 

\ 


Temperature 


Fic. The second moment (curve and the second central moment (curve CM) the 
line corresponding the first approximation. 

Fic. The second moment (curve and the second central moment (curve CM) the 
line corresponding the first approximation. 


SECOND APPROXIMATION FORMULAE FOR THE FIRST MOMENT 


The general expressions for the first and second moments the second 
approximation are rather complicated. have applied special cases only 
the general expression for the first moment (the general expression for the 
second moment the second approximation can found McMillan 
(1959)). 

One obtains 


a, 


8 


! 
| 
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Applying this expression the case the two absorption lines nickel 
fluosilicate, and using the assumptions formulated Section one obtains 
expressions which, when evaluated numerically, give the second approxima- 
tion curves Figs. and not reproduce here these expressions, 
because they are very long. They are the form 


4 


where depends the Boltzmann factors only, and depends the 
Boltzmann factors and the g-factors. 

should like emphasize the following features the curves given 
Figs. and 

The sign (Av), that the sign the shift the mean frequency the 
line relative ‘‘unperturbed”’ frequency determined for both lines 
the sign the exchange interaction constant Apart from this the tempera- 
ture dependence the mean frequency shift even qualitatively different 
for the two lines. This due the fact that the line corresponding 
originates the ground state single ion while the other line 
corresponds transition between excited states. very low temperatures 
the former line (transition from the ground state) which main interest. 
For this line the absolute value the mean frequency shift goes through 
maximum about the temperature decreases and becomes negligible 
below 0.4° 

The second approximation curves differ very little from the first approxima- 
tion curves. This means that one could obtain from the measured mean 
frequency shift very good estimate the exchange interaction constant 
since the first approximation the shift simply proportional 

One can prove quite generally that 


T>0 


the case transitions from the ground state any excited state for 
arbitrary value the effective spin assumption about the structure 
the crystal its shape need made nor necessary assume that 
commutes with the spin-Hamiltonian. 

The second central moment the case transitions from the 
ground state (as Fig. the first excited state can shown vanish for 
for arbitrary and quite generally. 

should add that have not seriously tried investigate the question 
what temperature the second approximation breaks down. course, 
possible that the qualitative features the temperature dependence the 
mean frequency shift remain unchanged higher approximations. 
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APPENDIX 
can arrive equation (3.6) the following way. 
Using definitions and notation explained Section one can easily see 
that the general expression for the /th moment the line shape function 


The summations over and are restricted those values for which 
> Ex. 

now evaluate this expression the defined 
Section others words, replace the energy values the spin 
system the Boltzmann factors the energy values the absence 
interactions: 


According the definitions given Section the matrix elements 


over only those values and for which This means that 
the numerator well the denominator equation (A.2) contains the 
factor After removing this common factor, one obtains 


| 

— 
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This expression becomes identical, except for the notation, with the right-hand 
side equation (3.6) all Boltzmann factors are replaced unity. sufficient 
condition for the validity such step course the validity the 
inequality 


which would, however, imply very high temperatures indeed. The reason why 
the temperature-independent equation (3.6) can used down tempera- 
ture the order 10° course due the fact that condition (A.4) 
sufficient without being necessary. rather obvious that sufficient for 
the temperature large compared with the differences between energy 
values single spin. formal proof this much weaker sufficient condition 
given the end Section see inequality (4.3). 


(A.4) (for all values and a’) 


DOMAIN STRUCTURE YTTRIUM IRON 


ARCHIBALD SMITH AND GERALD WILLIAMS 


ABSTRACT 


The magnetic domain structure thin slices yttrium iron garnet crystals 
has been studied transmitted polarized light. The transmitted light exhibits 
Faraday rotation and magnetic birefringence, and both effects have been used 
observe the domain structure. example structure samples having 
growth face one side presented for each the (110) and (112) planes. 
The structure for the (110) face similar that for ideal sample this 
plane. The structure for the (112) face differs from ideal sample due 
large strain anisotropy. 


INTRODUCTION 


has recently been found that many ferromagnetic oxides become optically 
transparent when the thickness reduced few thousandths centimeter. 
The transmitted light exhibits Faraday rotation and magnetic birefringence, 
and these effects have been used observe the magnetic domain structure 
(Dillon 1958; Sherwood, Remeika, and Williams 1959). contrast trans- 
parent films ferromagnetic metals, the oxide samples are sufficiently thick 
that the domain structure still characteristic the bulk material. this 
paper, some observations domain structure single crystals yttrium 
iron garnet are presented. These structures are analyzed briefly terms 
the appropriate magnetic theory. 

Yttrium iron garnet (abbreviated YIG) has previously been studied 
Dillon (1958). found that with samples ground both sides, the domain 
structure usually took the form irregular strips, which the magnetization 
was perpendicular the plane the sample, and antiparallel adjacent 
domains. This was interpreted being due strain the crystal caused 
the grinding and polishing operations. The strain results approximately 
uniaxial anisotropy perpendicular the surface which overrides the intrinsic 
cubic anisotropy. For convenience, this type structure will referred 
normal strain structure. Dillon’s work also suggested that the strain might 
less important only one side the sample were ground, the other being 
growth face. number samples were prepared investigate this point. 
Although some these exhibited the normal strain structure, the majority 
showed different types structure. Examples the two most interesting, 
occurring the (110) and (112) planes, are discussed below. The other 
structures observed were variations these two. 


EXPERIMENTAL 
2.1 Sample Preparation 
The crystals yttrium iron garnet were grown the method 
Nielsen and Dearborn (1958). melt and PbO the mole 
received March 25, 1960. 
Contribution from the Electronics Laboratory the Defence Research Telecommunications 
Establishment, Ottawa, Ontario. This work was carried out under Project No. D48-55-03-08. 
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ratio 0.525:0.440:0.035 was formed 1350° The temperature was then 
lowered C/hr about 950° Small crystals were obtained, usually 
not over along any direction, showing both (110) and (112) growth 
faces. The method grinding and polishing was similar that described 
Dillon and Earl (1959). The samples were mounted adjustable grinding 
jig, the piston which was left unclamped. With this arrangement, the 
sample was free ride over any slight irregularities the grinding surface 
without damage. Grinding was done with alumina glass, polishing with 
alumina phenolic-impregnated linen. Only one side the sample was 
ground, the other being growth face. The crystallographic plane the 
growth faces was determined X-ray techniques. 


2.2 Domain Observation 

When light which plane-polarized passes through ferromagnetic crystal 
parallel the magnetization, the plane polarization rotated. This 
type Faraday effect which arises from the internal magnetic field the 
crystal. general, the specific rotation varies rapidly with wavelength, but 
falls between and degrees/cm for YIG (Dillon 1959). example 
how the Faraday effect may used reveal domain structure, consider 
plane-polarized light passing through sample which has the normal strain 
structure discussed the Introduction. Because the alternation the 
direction the magnetization, the rotation will equal but opposite the 
two sets domains. the sample viewed through another polarizer 
right angles the first, the two sets domains will have equal intensity. 
The domain walls will appear black, however, because much their magnetiza- 
tion perpendicular the line sight and hence does not cause any rotation. 
Contrast between the two sets domains can obtained uncrossing the 
two polarizers slightly. YIG also exhibits magnetic circular dichroism (Dillon 
1959): the absorption right- and left-handed circularly polarized light 
slightly different, and the sign the difference depends whether the light 
travelling parallel, antiparallel, the direction the magnetization. 
This effect can also used detect domain structure but does not give 
good contrast the first method. 

When plane-polarized light passes through magnetic crystal perpendicular 
the magnetization, general emerges elliptically polarized due the 
magnetic birefringence. The ellipticity varies from zero when the magnetization 
either parallel perpendicular the plane polarization, maximum 
when the magnetization 45° the plane. The magnetic birefringence 
very weak, giving retardations only few millimicrons for transparent 
samples. This effect provides means determining the direction the 
magnetization the plane perpendicular the line sight, with ambiguity 
90°. 

The apparatus for domain observation consisted standard polarizing 
microscope* with photographic attachment, and two small electromagnets 
which could fixed the rotating stage. One magnet gave field perpendic- 


*Model LI-2, Bausch and Lomb Optical Co., Rochester, N.Y. 
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ular, the other parallel, the line sight. When using these magnets, the 
regular condensing lens was replaced suitable strain-free objective lens. 
Fields were measured with Hall-effect gaussmeter. The retardation elliptic- 
ally polarized light was measured the quarter-wave plate method (Goranson 
and Adams 1933). order carry this out, quarter-wave plate was mounted 
below the stage specially-constructed holder, and the lower polarizer was 
provided with finely-graduated angular scale. The latter was also used 
measure the Faraday rotation. 


RESULTS AND DISCUSSION 


3.1 Domain Structure the (110) Plane 

The domain pattern Fig. was characteristic strain-free samples cut 
from (110) growth faces. less extensive example this pattern was shown 
Dillon (1958). The sample exhibited only magnetic birefringence and thus 
the magnetization lay the plane the sample. When the sample was 
viewed between crossed polarizers only, the contrast between the domains 
was very low. Therefore, taking the photographs this sample, the contrast 
was increased using quarter-wave plate. Normal strain structure occurred 
along two edges the sample. The vertical and horizontal crosslines Fig. 
are parallel the and directions, respectively, and the same time 
indicate the position the crossed polarizers. The crystallographic directions 
are explicitly marked Fig. which idealized drawing Fig. 

Two types domain boundaries could discerned this sample. The first, 
which will called type was clearly visible the boundaries between the 
light and dark areas, and lay the [001] and [110] directions. The second, 
type took the form fine white black lines cutting through either light 
dark region, and ending the right-angle junctions type boundaries. 
The acute angle between type boundaries and the [001] direction was 
approximately 57°. several cases, these boundaries were rather broad, and 
appeared inclusions the normal strain structure. The type boundaries 
were barely visible the original photograph Fig. la, and may lost 
reproduction. For clarity, the boundaries have been shown schematically 
Fig. 

The change intensity which occurs crossing type boundaries indicates 
that these are neither 90° nor 180° walls. Conversely, type are either 90° 
180° walls. clear from Fig. that the only consistent interpretation 
type 114° and 66° walls the [001] and [110] directions, respectively, 
and let type 180° walls. order prevent the formation poles 
within the sample, the 180° walls must terminate junctions 114° and 66° 
walls, i.e. right-angle corners. will seen that the four orientations 
the magnetization, and the 180° walls, are nearly parallel the four (111) 
directions the plane. 

When field was applied along one the (111) directions the plane, one 
set domains, presumably that with the magnetization the field direction, 
grew the expense the others. The type boundaries moved right 
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angles their length, such way maintain 66°, 114°, and 180° walls 
until the domains were swept out. This occurred field about oersteds. 
However, few hazy, wedge-shaped domains were left the sample this 
field, and about oersteds was required remove all traces structure. 

The behavior with the field the direction was more complex. 
Initially, two sets domains, one light and one dark, began increase 
size, with the boundaries moving before. These two sets formed strips across 
the sample perpendicular the field, shown Fig. 1b, that only 114° 
walls remained. the field was further increased, the strips became narrower 
and their number increased. the same time, another type domain with 
intermediate intensity began appear, again the form strips perpendic- 
ular the field. These domains can seen Fig. still higher fields, 
the contrast between the light and dark strips decreased, and finally they 
merged with the intermediate domains field 100 oersteds. However, 
some traces structure remained until field 140 oersteds was reached. 

With the field the direction, the behavior was similar the 
previous case. Strips again formed perpendicular the field, but the inter- 
mediate domains entered much earlier this case, and clear-cut strip pattern 
like that Fig. did not occur. rather hazy pattern with wedges instead 
strips was obtained, somewhat similar Fig. The sample became nearly 
uniform appearance field oersteds, although some traces struc- 
ture remained until the field reached oersteds. 

The behavior the two latter cases may understood taking the 
magnetization the light and dark strips along the two (111) directions 
which were closest the field direction. order saturate the sample, the 
magnetization had rotated until was parallel the field. This was 
evidently occurring two ways: gradual rotation within the black and 
white strips, and sudden rotation, leading the formation the inter- 
mediate strips. 

The domain structure predicted for ideal sample the (110) plane will 
now briefly described. For further details and references, the review article 
Kittel and Galt (1956) may consulted. YIG has negative cubic aniso- 
tropy, and the four (111) directions the (110) plane are directions easy 
magnetization. These are separated angles 70°32’ and 109°28’. The 
simplest domain structure expected occur when the sample equilateral 
parallelogram with the sides the (111) directions. consists four quad- 
rants, each with the magnetization one the easy directions, such way 
give flux closure within the sample. The boundaries are the [001] and 
directions, ending opposite corners the sample and crossing right 
angles. When field applied along (111) direction, the boundaries are 
expected move perpendicularly their length until the three unfavorable 
domains are eliminated. When field applied along either the [110] 
direction, the structure changes series strips perpendicular the field, 
with the magnetization alternating the two favorable easy directions. The 
number strips increases with increasing field, and the magnetization rotates 
towards the field. The fields required for saturation are approximately 
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and the [001], and [110] directions, 
respectively, where the magnetization, the anisotropy constant, and 
the demagnetization factor the plane the sample. 

The observed structure and its behavior under applied field are con- 
sistent with these predictions. The differences can readily understood 
the basis the irregularities the crystal structure and shape the sample. 
addition, the occurrence the intermediate domains when the field 
the [001] [110] direction probably connected with the low value the 
demagnetization factor The value estimated from the dimensions 
the sample only 0.013. The predicted saturation fields are 22, 106, and 
oersteds the [001], and [110] directions, respectively, with 
1750 K/M —42 oersteds, and 0.013. These are reasonable agree- 
ment with the experimental fields 25, 100, and oersteds required remove 
the gross domain structure. 


3.2 Domain Structure the (112) Plane 

The sample Fig. was cut from (112) growth face. displays two 
domain patterns: narrow and. broad strips, both parallel the [111] direction. 
(The circle near the center this figure bubble the mounting wax.) 
before, the cross hairs indicate the position the polarizers. With the 
polarizers exactly crossed, the contrast between the narrow strips was 
maximum when they were oriented 45° the polarizers, and was practi- 
cally zero when they were parallel the polarizers. Under the same conditions, 
the contrast between the broad strips was very low, regardless orientation. 
was found that better contrast could obtained uncrossing the polari- 
zers few degrees, and this has been done Fig. The region with the 
narrow strips exhibited permanent birefringence, independent the magneti- 
zation. The direction vibration the slow wave was parallel the [111] 
direction. The retardation was millimicrons, about five times the expected 
magnetic birefringence. The region with the broad strips exhibited weaker 
permanent birefringence, being about the same the magnetic birefringence. 
Closure domains could just discerned this region, near the edge the 
sample. They took the form spikes protruding into the strips. The thickness 
the sample this region was about half that the region the narrow 
strips. 

The behavior this sample with external field applied was investigated 
for three orientations the field. (1) Field the plane, parallel the strips. 
The domain walls moved perpendicularly themselves, that one set 
strips each region expanded, while the other shrank. The narrow strips 
merged field 550 oersteds, the broad strips oersteds. (2) Field 
the plane, perpendicular the strips. Initially the region with the narrow 
strips brightened slightly, and then the narrow strips formed zigzag pattern 
shown Fig. 2c. the field was increased, the contrast between the strips 
decreased. the same time the angle between the strips and the field direction 
decreased, and became almost zero when the pattern disappeared 600 
oersteds. The widths the light and dark strips appeared remain equal, 
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but decreased somewhat the field increased. The broad strips were unstable 
with this field orientation, and again disappeared relatively low field. 
(3) Field perpendicular the plane. the domain walls moved 
case (1). 600 oersteds the broad strips became unstable, and irregular 
set domains with the magnetization perpendicular the plane began 
appear. higher fields, the broad strips were engulfed these irregular 
domains, which grew until they merged together 1750 oersteds. The narrow 
strips merged 930 oersteds. the field was further increased, the Faraday 
rotation this region increased until reached maximum value 2300 
oersteds. This indicated that the magnetization was being rotated into 
position perpendicular the plane. 

With negative cubic anisotropy, there one easy direction the (112) 
plane, namely [111]. Ideally the zero-field domain structure would expected 
consist strips, with the magnetization the plane the sample and 
antiparallel adjacent strips. The results for case (1) above are consistent 
with such structure. However, neither the Faraday rotation nor the magnetic 
birefringence would provide any contrast between the domains. The existence 
such contrast shows that the magnetization has component perpendic- 
ular the plane. This consistent with the results for case (3), which also 
require perpendicular component, antiparallel adjacent domains. Thus the 
magnetization inclined the plane, and this, together with the existence 
the permanent birefringence, requires that the crystal strained, with 
the strain axis inclined the plane. The strain must induce approximately 
uniaxial anisotropy which overrides the cubic anisotropy the region the 
narrow strips, but much weaker the region the broad strips. explain 
the existence strips the direction, the projection the strain axis 
the plane must lie this direction, which consistent with the orientation 
the optical axes. 

The contrast between the narrow strips can explained terms the 
circular dichroism which occurs YIG. Consider first the case which the 
strips are 45° the polarizers. When linearly polarized light passes through 
the sample, emerges with small circular component due the weak bi- 
refringence. The linear component may broken down into two equal circular 
components rotating opposite directions. The small circular component can 
added one these, giving two circular components with slightly different 
amplitudes. the magnetization has component parallel the light 
path, one the circular components will absorbed more than the other, 
due the circular dichroism. For one direction the difference between 
the circular components will reduced, and for the opposite direction, 
increased. The resulting difference intensity for the two directions 
detected the second polarizer. When the strips are parallel the polari- 
zers, the light emerges linearly polarized, and the absorption affects all domains 
equally. For the broad strips, the birefringence apparently too weak 
provide appreciable contrast this mechanism. However, this can over- 
come making use the Faraday rotation associated with i.e. 
uncrossing the polarizers slightly, Fig. 2a. 
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The fact that the magnetization was inclined suggested that the sample 
tilted under the microscope. apparent depth measurement showed that the 
index refraction YIG was about 2.2. Therefore, order increase the 
internal ray angle, the sample was immersed methylene iodide, which has 
index 1.75. The methylene iodide was placed small container with 
windows perpendicular the line sight. The experiment was carried out 
with the tilt axis perpendicular the strips, which were 45° the polarizers. 
was expected, the contrast between the strips varied with the angle tilt. 
Positions minimum contrast existed for both sets domains. With the 
sample suitably tilted, another set broad strips becomes visible, super- 
imposed the narrow strips. This illustrated Fig. 2b. noting the 
optimum focussing positions, was established that the two sets strips 
existed separate layers, with the narrow strips closer the growth face. 
can seen Fig. 26, smooth transition occurred between the two sets 
broad strips. 

The position minimum contrast between the narrow strips occurred 
tilt 20°. The contrast could restored uncrossing the polarizers slightly, 
indicating that the internal ray direction was parallel the magnetization. 
For the first set broad strips, this position occurred 4°, but the contrast 
was not affected uncrossing the polarizers. This indicated that the internal 
ray direction was perpendicular the magnetization. For the second set 
broad strips (under the narrow ones), the corresponding position occurred 
0°. The inclination the magnetization the plane was calculated 
75°, 3°, and 0°, respectively, for the three cases. 

may concluded that the highly strained material occurs layer 
base unstrained material. The region without the narrow strips essen- 
tially continuation the base layer, but weak strain has been introduced. 

The behavior with applied field can easily least qualita- 


Fic. schematic representation Fig. The solid and dotted lines represent type 
and boundaries, respectively, and the arrows indicate the direction the magnetization. 
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tively, except when the field perpendicular the strips (case 2). The initial 
brightening the sample can attributed the saturation the unstrained 
base layer. The reduction contrast between the narrow strips indicates that 
magnetization the strained layer rotating into the plane. It-appears that 
the magnetization remaining parallel the domain walls. order get 
net component magnetization the field direction with this arrangement, 
the component magnetization the plane must the same direction 
for all parallel strips. However, this can only achieved the magnetization 
every second strip has flipped over position with the opposite tilt. The 
formation the zigzag pattern permits the boundaries remain simple 
variant the original 180° walls, instead becoming more complex walls. 
Thus this analysis requires that the increase anisotropy energy due the 
flipping the magnetization more than offset the avoidance complex 
domain walls. 

The permanent birefringence was found all the four samples which 
were cut from (112) growth faces. addition the present sample, one 
other was completely analyzed. There were layers this case, and the 
magnetization was inclined angle 60° throughout the sample. The 
direction vibration the slow wave was again parallel the direction. 
permanent birefringence was also found one sample cut from (110) 
growth face. clear that the strain giving rise the birefringence was 
introduced during the growth the crystals. Although the present evidence 
rather limited, appears probable that all (112) growth faces are strained. 
This would consistent with the observation Nielsen and Dearborn (1958) 
that (112) faces predominate when the crystals are grown rapidly. 


CONCLUSIONS 

has been found that samples having growth face one side are less 
likely develop strain from the grinding and polishing operations. The domain 
structure samples which are free from this type strain depends the 
particular crystallographic plane. Those cut from (110) faces have structure 
which similar that expected for ideal sample this plane. Those cut 
from (112) faces have unusual structure, due large strain anisotropy which 
introduced during the growth the crystals. 


wish thank Mr. Peter Buratynski for grinding and polishing the 
samples, and Mr. Jack MacDonald for making some preliminary observa- 
tions. 
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THE POLAR RADIO BLACK-OUT THE IONOSPHERE 
JULY, 


CHAPMAN 


ABSTRACT 


sounding station data, particularly were used study the 
ionospheric storm July, 1958, which shown have consisted clearly 
two stages. Polar cap absorption began shortly after importance flare 
was observed 0039 U.T., and magnetic storm and auroral absorption event 
followed 0749 U.T. July 

The polar cap absorption spread towards the geomagnetic equator, but the 
use higher than first order magnetic dipole terms was found necessary order 
arrange the data properly latitude. The absorption began some hours earlier 
northern stations than most southern stations, the main cause the dif- 
time onset being thought due the presence absence 
sunlight. 

The inference drawn that ionospheric absorption measured about Mc/s 
caused ionization 50- 80-km heights, which case the data suggest 
that electrons first precipitated the earth above 83° corrected geomagnetic 
latitude, and then protons and helium ions lower latitudes. The cutoff latitude 
was 59°, corresponding particle magnetic rigidity ev/c. This 
storm was compared with the more severe event February, 1956, and showed 


INTRODUCTION 


The phenomenon prolonged periods severe radio absorption the 
polar regions has been known for some years (Obayashi 1959; Leinbach and 
Reid 1959; and Bailey 1959). Until the however, there were insufficient 
data available for adequate description. now clear that during major 
ionosphere storm this type, there absorbing blanket over the whole 
the polar cap, which has led the name cap being 
applied. Such event occurred during the period July July, 1958, 
coinciding with Special World Interval (July (Harang and Troim 
1959; Hultqvist, Aarons, and Ortner 1959; and Hultqvist 

The event presumed have begun shortly after midnight July, 
1958, with importance solar flare which was observed from 0039 
0324 U.T. the Tokyo Astronomical Observatory Mitaka and from 0040 
0214 U.T. Honolulu, Hawaii, and which reached its maximum intensity 
about 0103. The identification the time beginning this storm 
reasonably certain, the only possible confusion arising from importance 
flare, observed Mitaka from 0032 0117 U.T. July. Both these 
flares were situated near the central median the sun about degrees 
north latitude. There were flares importance greater than observed 
during the previous day, which was international magnetically quiet 
day for the month. 

Major radio noise bursts were recorded shown the following table 
(CRPL Series, and High Altitude Observatory, Boulder). 
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TABLE 
Solar radio noise emission 


Time duration Frequency Location 
(Mc/s) observatory Remarks 

10,000 Ottawa, Canada 
0029-0058 458 Boulder, Colo. 
0029-0131 167 Boulder, Colo. 
Fort Davis, Texas Type III fast drift burst 
0028-0029 100-580 Fort Davis, Texas Continuum, intensity 
0052-0147 100-580 Fort Davis, Texas Type IV, continuum, intensity 
0055-0152 10,000 Ottawa, Canada 
0058-0207 458 Boulder, Colo. 
0131-0214 167 Boulder, Colo. 


These data appear fall into two groups, one beginning about 0027 U.T., 
and the second 0052 U.T. The second group might associated with the 
flare 0039 U.T., and the first with the flare observed earlier. The 
maximum intensity the noise bursts 0104 coincided with the peak 
solar flare activity the later flare. 

These flares and noise storms were also associated with short-wave fade-outs 
(SWF), particular widespread slow SWF major importance from 0025 
0230 hours U.T. 

The first evidence the onset the storm was gradual increase iono- 
spheric absorption observed the ionospheric stations within the northern 
auroral zone, north 65° geomagnetic latitude. 0500 U.T., severe 
black-out was effect over the whole the north polar cap, followed 1100 
U.T. similar event the southern polar region. This black-out continued 
with letup until the beginning July. 

The polar cap black-out was accompanied relatively minor magnetic 
activity. 0749 U.T. July severe world-wide geomagnetic storm 
began with sudden commencement (SC) (Shapley 1959). Although the north 
polar regions are largely sunlit this time the year, the evidence suggests 
that severe auroral storm followed. 

Cosmic ray monitors the northern hemisphere recorded Forbush 
decrease, beginning between 0900 and 1000 U.T. July and continuing 
until July 12. 

The purpose this paper present the data this particular polar 
radio black-out, and describe some significant features, the hope that 
similar features can recognized other storms this type. 


IONOSPHERIC SOUNDING DATA 


The data used this study were obtained the ionospheric sounding 
network established during the I.G.Y., from the stations shown Fig. for 
the northern hemisphere, and Fig. for the southern hemisphere. 

The principal data used this study were f-plots, showing the frequencies 
reflected from the ionosphere function time day. The progress the 
storm can seen Fig. which shows the f-plots for six Canadian ionospheric 
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CAMPBEL( ISLAND 


© SOUTH GEOMAGMETIC POLE 


Fic. stations the southern hemisphere. 


stations, for the period July July 12, arranged order geomagnetic 
latitude. The ordinate shows the critical reflection frequency megacycles 
measured every minutes, the abscissa the time day. 

The beginning the storm can seen Resolute Bay, Baker Lake, 
Yellowknife, and Churchill gradual increase the minimum frequency 
echo from the ionosphere denoted vertical line extending from 
the base line and terminated dot. solid vertical line extending the 
top the graph indicates complete black-out, with echoes being received 
from the ionosphere any frequency. should noted particularly that the 


f-plot from Ottawa showed unusual absorption following the time the 


flare, until the sudden commencement hours later. 

The extent the polar black-out can seen Figs. and showing the 
north and south polar regions respectively with data from the various stations 
arranged order geomagnetic latitude. apparent that the storm began 
later the south than the north, except for short period black-out 
Wilkes Station 0500 U.T. July which corresponds with the time 
black-out northern stations comparable geomagnetic latitude. 
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MAGNETIC, AURORAL RADAR, AND COSMIC RAY DATA JULY 


order show the two main phases the storm, magnetic, auroral radar, 
and cosmic ray data are shown Fig. for the period from July July 10. 

The main features the variations the horizontal component the 
magnetic field shown for Agincourt were similar Baker Lake, Resolute Bay, 
and Alert. The fluctuations the magnetic field prior 0749 July were 
relatively minor. The severe magnetic storm which followed the sudden com- 
mencement caused large and rapid magnetic field fluctuations for about 
hours, which were similar the four magnetic observatories. 

The v.h.f. auroral radars Ottawa and Resolute Bay showed returns 
prior the sudden commencement. However, strong auroral echoes were 
observed north the Ottawa station the time the sudden commencement 
and continued off and until July 10. The geometry reflection such that 
these echoes come from the vicinity the auroral zone maximum. Resolute 
Bay, radar aurora echoes were observed except for weak echo present for 


The cosmic ray data show typical Forbush decrease beginning about the 
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same time the sudden commencement. The cosmic ray count did not return 


lar, normal until about July 12. should noted that there was noticeable 
10. increase cosmic ray intensity before the Forbush decrease, there was 

the the time the February 23, 1956, event. The peak 0800 Resolute Bay 

ay, not considered significant this regard. 

ere The storm may therefore divided into two distinct phases. The first began 

with the solar flare activity, the second with the sudden commencement the 
geomagnetic storm. Because these two phases are separated period 


hours, possible study each one separately. The intent this paper 
examine detail only the first phase. 


TIME ONSET POLAR BLACK-OUT 

hat The results detailed study the time onset polar black-out are 
lute shown Fig. which gives the deviations from the monthly medians 
for similar manner that used Shapley and Knecht (1958) for the 
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February, 1956, event. The times last normal sounding, first possible 
and first definite increase and the first black-out are shown Table 
for later analysis. 

These data show the gradual increase absorption from the time the 
flare, until complete black-out about 0500, northern stations. Evidence 
the short-wave fade-out (SWF) can seen some the stations sun- 
light increase immediately after the flare, before the polar cap 
absorption began. The initial black-out extended St. John’s, 58° geo- 
magnetic latitude. Three stations, however, Table II, Inverness, Kjeller, 
and Upsala are clearly out order, since black-out did not occur these 
stations until after the sudden commencement July 

the south the black-out extended north Ellsworth (67°), and there 
was evidence increased absorption Campbell Island (57°) during daylight 
hours. The late appearance absorption Ellsworth is, however, difficult 
explain for station this high geomagnetic latitude. 

The solution these inconsistencies connected with the latitude the 
stations found use more precise values for the geomagnetic lati- 
tude. Studies have been made Hultqvist (1958, and Quenby and 
Webber (1959) using fifth and sixth order harmonics the geomagnetic field 
order obtain better representation the effective magnetic latitude. 
These authors have shown that application such corrections lead more 
accurate representation auroral isochasms. 

The calculations Hultqvist (1958) have been used compute corrected 
geomagnetic latitude tabulated Table II. 

The inconsistencies mentioned above disappear when the stations are 
considered order corrected geomagnetic latitude. Inverness, Kjeller, and 
Upsala are placed south St. John’s, which 59° corrected latitude the 
southern limit the polar black-out. 

the southern hemisphere, Ellsworth (62°) and Campbell Island (61°) 
compare well with Winnipeg (61°) and St. John’s (59°) when times black-out 
are compared, taking into account the amount daylight each station. 
The cutoff magnetic latitude for this storm, therefore, about 59°. 

was shown Reid and Collins (1959) that polar cap absorption greater 
during daylight hours. presumed that sunlight causes the dissociation 
negative ions, maintaining high electron densities the lower ionosphere. 
This effect can seen Winnipeg, St. John’s, Wilkes Station, and Campbell 
Island. The early onset black-out Wilkes Station presumably partially 
due the short period daylight that station, which, though high geo- 
magnetically, 66.4° geographic latitude. The second onset black-out 
Wilkes about 1100 hours roughly corresponds with the time black-out 
Pole Station, America, Cape Hallett, and Byrd Station the south, 
and Lycksele, Meanook, Winnipeg, and St. John’s the north. 

therefore suggested that the main reason for the late onset polar cap 
absorption southern stations due these stations being darkness. The 
alternative explanation, that the inclination the earth’s magnetic field dipole 
axis away from the sun this time year reduced the incoming particle flux 
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seems less likely, since attempt correlate local time particle precipita- 
tion with opening the Stérmer cone, using the graphs produced Alpher 
(1950), was inconclusive. 

The onset black-out within hour the sudden commencement, 
Winnipeg, St. John’s, Ottawa, and Campbell Island presumably shows either 
the beginning auroral absorption stations outside the auroral zones, 
the equatorial extension polar cap absorption due the distortion the 
earth’s magnetic field. 


LATITUDE VARIATIONS TIME FIRST ABSORPTION 


The times first increase absorption the stations listed Table 
were rearranged according corrected geomagnetic latitude, and plotted 
Fig. principle, the time first increase should relatively more 
reliable than, for example, time first black-out, since the latter function 
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Fic. first increase arranged order corrected geomagnetic latitude, 
July, 1958. 


equipment sensitivity, antenna pattern, ionospheric critical frequency, and 
on. The increase relative measurement based average character- 
istics the equipment and the ionosphere each frequency. 

The early beginning absorption highest latitude stations and the 
subsequent spread absorption lower latitudes are clearly evident. the 
highest latitude stations, the measurement the precise time onset 
absorption made somewhat more difficult the short-wave fade-out, which 
evident number the stations sunlight. should noticed particu- 
larly that the first possible increase absorption Wilkes Station Antarctica 
occurred approximately the same time the increase northern 
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TABLE 
Geo- Corrected First 
dipole magnetic normal increase increase 
Station Abbr. Operated by: latitude 
Thule USS. Signal Corps 84°N 0130 0200 0230 
Eureka Canada D.R.B. (0000) (0100) (0200) 
Alert D.R.B. 85°N. 87°N. 0145 0200 0200 
Resolute Bay Canada D.O.T. 85°N 0100 0215 0215 
Clyde River D.R.B. 82°N. 82°N. 0245 0300 0300 
Fletchers Signal Corps 83°N. 0200 0215 0215 
Research Est. 
Churchill D.O.T. 68°N. 0645 0700 0715 
Barrow 69°N. 0300 0345 0345 
Yellowknife Canada D.R.B. 0300 0315 0315 
Tromso Norwegian Def. (0300) (0400) (0400) 
Research Est. 
Kiruna Inst. 65°N. 0300 0400 0400 
Defence, Sweden 
Lycksele Inst. Nat. (0900) (1000) (1000) 
Defence, Sweden 
Meanook Canada D.R.B. 62°N. 62°N 0915 0930 0930 
Inverness U.K. 61°N. 57°N 
Research Est. 
Winnipeg D.O.T. 61°N. 1000 1030 1030 
Defence, Sweden 
St. John’s U.S. Signal Corps 1200 1215 1215 
Ottawa Canada D.O.T. 59°N 
phys. Obs. 
phys. Obs. 
Ellsworth 1800 1915 2100 
Little America 0445 0500 0500 
Cape Hallett Geo- (0500) (0700) (0700) 
phys S 
Wilkes Station 0300 0315 0415 
Pole Station 0530 0615 0645 


stations the same latitude. The delay time first absorption dark 
stations was the order hours except Churchill which showed some 
absorption and then recovered, which made assessment the time first 
increase very uncertain. 

striking observation this graph that one group data appear 
clustered above 83° effective magnetic latitude, about 0200 U.T. Absorption 
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Duration July 
Black-out range 

0500 6.0-7.5 
0515 5.2-6.3 
0500 116 4.3-5.9 
0515 5.0-6.8 
0500 5.2-6.9 
0415 5.1-6.7 
0500 4.6-7.0 
0545 5.0-7.3 
0430 3.9-7.2 
0815 5.4-7.0 
0345 5.0-6.4 
0545 5.0-6.5 
0500 
0700 4.5-7.3 
0515 5.2-6.6 
1400 4.5-7.7 
1030 4.5-6.7 
5.5-8.2 
1100 3.8-6.8 
1215 4.6-7.6 
3.8-7.3 
4.6-8.0 

3.1-13 

2.8-13 
2100 1.9-5.2 
1045 2.4-9.5 
1045 2.2-8.0 
1300 2.9-8.4 
0445 62+ 3.3-8.0 
1000 62+ 3.6-8.2 
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Absorption during 
July 


Black-out progress 
Black-out progress 
Black-out progress 
Black-out progress 
Black-out progress 
Black-out progress 
Black-out progress 


Black-out progress 
Black-out progress 
Black-out progress 
Black-out progress 
Black-out progress 
Black-out progress 


Black-out progress 


Black-out progress 
Black-out progress 


Black-out progress 

Black-out began 1500 
July 

Black-out began 1500 
July 

Absorption high, 
black-out intermittent 

Black-out began 1500 
July 

Intermittent until 
black-out 0815 July 

Black-out began 0900 
July 

Intermittent black-out 
began 1415 July 

Normal 


Black-out began 0900 
July8 

Black-out progress 

Black-out progress 

Black-out progress 

Black-out progress 


Intermittent black-out 
Blackout progress 


Remarks 


f-Plot, hourly values only 


f-plot, hourly values only 


Equipment unserviceable 
f-plot, hourly values only 


Equipment unserviceable 

f-plot, hourly values only 

values only 


f-plot, hourly values only 
9 hr 


f-plot, hourly values only 


f-plot, hourly values only 


f-plot, hourly values only 


did not occur any station lower latitude before 0300 U.T. The lowest 


latitude which absorption was evident before the sudden commencement 
was 59°, pointed out previously. 

For shows the time onset black-out. Again the 
blacked about hours before stations darkness, except 


for 


low 62°. 


0200) 
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0300 
0300) 
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Fic. Time onset black-out, arranged order corrected geomagnetic latitude, 
July, 1958. 
MECHANISM ABSORPTION 


increase electron density the lower ionosphere (50-80 km) 
undoubtedly the cause the increased absorption. The coefficient absorption 
given the following equation: 


where and are the electronic charge, the electron density, the collision 
frequency, and the radio wave angular frequency. The usual approximations 
have been made neglecting the magnetic field and assuming that the index 
refraction not significantly different from unity. 

constant, the absorption coefficient reaches maximum the 
height which decreasing heights the collision frequency becomes 
much greater than and that absorption soon becomes 
unimportant. The variation collision frequency with height was measured 
rocket from over Churchill during polar black-out Kane 
(1959). The measured values, extrapolated down and 
are shown Table III. 

From this table evident that absorption about Mc/s 
will small below about 10w) and will reach its maximum about 
w). Measurements Seddon and Jackson (Friedman 1959) 
the same rocket flights above Churchill, when, however, the electron density 
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TABLE 


Collision frequencies 


was not constant, showed that absorption 7.75 Mc/s began about 
the rocket’s upward flight, reached its maximum km, and 
decreased rapidly above this height. 

can inferred from this argument that ionospheric sounding method 
observation absorption selective. The measurements reported this 
paper emphasize the effects which produce ionization heights from approxi- 


SOLAR PARTICLE PRECIPITATION 


The causative agent for the ionization these heights the ionosphere 
believed the precipitation charged particles from the sun produced 
some mechanism, yet unspecified, beginning with flare major 
importance (Leinbach and Reid 1959; Hultqvist Reid and Collins 
1959; and Obayashi 1959). 

The particles lose their energy collisions the atmosphere. Each particle 
produces very large number ion pairs ion pairs for 10-Mev 
proton), the majority being produced the last few kilometers air where 
the solar particle losing energy rapidly. 

The particles are constrained move Stérmer orbits the earth’s 
magnetic field since the particle flux extremely low. The lack reaction 
the earth’s magnetic field during the first hours polar cap absorption 
event confirms this low flux density. consequence reasonable use 
the undisturbed magnetic field values compute the magnetic cutoff rigidity 

The magnetic cutoff rigidity, corresponding magnetic latitude given 
(1955) 


where the magnetic dipole moment the earth, and the radius 
the earth. The values for the appropriate latitude are plotted the 
right-hand ordinate axes Fig. Also plotted are the kinetic energies and 
penetration depths the atmosphere for particles corresponding rigidity. 
Proton and helium ranges from the summary Whaling (1958), electron 
ranges after Heitler (1954), and atmospheric parameters after Nicolet (1959) 
were used. Following Bailey (1959) has been assumed that the helium ions 
are singly charged coming from the sun, but that they lose their second 
electron through collision entering the atmosphere, and penetrate depths 
appropriate for alpha-particles. 

apparent that the 83° cutoff Fig. could interpreted due 
electrons with energies about Mev, providing the earth’s magnetic field 
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essentially dipole character these latitudes. difficult see how 
heavier particles could responsible for this effect, since proton He+ 
ion flux which would penetrate into the atmosphere, would also 
evident much lower geomagnetic latitudes. The solar radio noise burst 
following the solar flare has been widely interpreted due radiation from 
energetic electrons, and their precipitation into the earth’s atmosphere 
not entirely unexpected. 

There some indication Fig. break 70°, which might due 
protons, while the low latitude cutoff 59° suggests that present 
the particle flux. 

well known that protons are present during polar cap absorption 
events (Leinbach and Reid 1959). The presence has been discussed 
Hultqvist (1958) and balloon observations cosmic particles July 
have been interpreted possible evidence solar alpha-particle flux Meyer 
(1959). is, course, well known that about 15% the galactic cosmic 
ray flux due ions. 


COMPARISON WITH FEBRUARY, 1956, EVENT 


interesting compare the data from July, 1958, with that from 
the similar study Shapley and Knecht (1958) the February, 1956, 
event. These authors analyzed data from ionospheric stations the northern 
hemisphere. 

flare 0334 U.T. February presumed have initiated the 
earlier event, which, however, occurred when the north polar region was 
darkness. The data Shapley and Knecht have been replotted Fig. 

The delay between the initiating flare and the sudden commencement was 
hours, and the duration black-out was about hours for the 1956 
event. The low latitude limits for the ionospheric polar cap black-out 
July and February events were 59° corrected magnetic latitude (St. John’s) 
and 58° (Kjeller) respectively, corresponding magnetic rigidities 9.6 
and electron volts. However, during the daytime when polar cap 
absorption more efficient, the effect February extended 54° (Fort 
Monmouth) indicating that particles with rigidities electron volts, 
were present. Protons helium ions this rigidity would penetrate very 
low altitudes, and either the flux was very large, heavier ions were respon- 
sible for this absorption. Cosmic ray data from equatorial latitudes showed, 
course, the presence particles much greater energy. 

The more rapid extension polar absorption lower latitudes, and the 
very marked increase cosmic ray intensity, show that the February 
event was the more energetic the two, even though the duration black-out 
was less. Otherwise, there are marked similarities between these two events 
which show that they were both the same class phenomenon. 


CONCLUSION 


From this study the July 1958, event, clear that polar cap absorp- 
tion caused particle precipitation into the earth’s atmosphere which 
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Fic. 10. Time first increase arranged order corrected geomagnetic latitude, 
February, 1956 (after Shapley and Knecht). 


produced increased ionization and absorption the lower ionosphere 
heights the order The earth’s magnetic field modified the 
selected flux particles, that measurable flux did not penetrate 
magnetic latitudes below 59°, corresponding particle rigidity 
ev/c. The necessity for use better approximation than the dipole field 
determine the precipitation trajectories particles has been shown. 

The later onset the disturbance the southern hemisphere apparently 
was due the absence sunlight the south polar regions, rather than 
peculiarities the particle orbits. The data suggest that electrons first preci- 
pitated into the ionosphere causing absorption latitudes higher than 83° 
corrected magnetic latitude, and that protons and helium ions then precipitated 
lower latitudes. 

The polar cap absorption event was followed sudden magnetic com- 
mencement auroral absorption, and Forbush decrease. The effects 
were very similar those the February, 1956, event but were longer 
duration, particularly the northern hemisphere, though the data suggest 
that the particles were not quite energetic the earlier occasion. 
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VELOCITY AUTOCORRELATIONS CHARGED PARTICLES 
MAGNETOIONIC MEDIUM WITH APPLICATIONS 
TURBULENT 


ABSTRACT 


system charged particles slightly ionized medium considered 
subject (1) collisions with members the neutral species, (2) constant 
external magnetic field, and (3) fluctuating force field, either external 

representing the collision forces. the assumption that their motions are 
satisfactorily described the Langevin equation, the cross-correlation functions 
time between velocity components these charged particles are calculated. 
These functions may used, described elsewhere the author, describe 
the self-diffusion the charged particles. The cases treated are: purely random 
external forcing, forcing exponentially correlated (Markovian) forces, and 
forcing random series pulses corresponding collision forces. 


INTRODUCTION 


The term used here slightly extended sense: 
what meant diffusion which the motion the diffusate particles 
coherent over finite space and time intervals. This contrasts with the com- 
pletely random motion implicit the classical (Fickian) diffusion equation. 
the particular case which the diffusion due turbulent flow the 
medium containing the diffusate particles, such the diffusion dissolved 
salts the sea, the diffusion turbulent the usual sense. All diffusion ex- 
cepting the limiting classical case turbulent. the present article, the 
time correlation functions the velocity components will calculated for 
few simple cases. The spatial correlations will neglected here because (1) 
not known just how they influence the quantitative description the 
diffusive process, (2) not clear how they may calculated for the physical 
systems considered here, and (3) the systems treated here are such that the 
effect the spatial correlations assumed negligible. The characterization 
the turbulence will taken the velocity autocorrelations first 
order, i.e., quadratic the velocities, although recognized that complete 
characterization must include correlations all orders. These higher order 
correlations would probably have only small effect the diffusive process 
and, any case, even the first order correlations have not yet been taken 
properly into account. 

well-known kinematical theorem the mean square spread the diffu- 
sate originating point depends upon the velocity autocorrelation function 
(VAF) according 


(1) 


received March 14, 1960. 
Contribution from the Hughes Research Laboratories, Hughes Aircraft Company, Malibu, 
alifornia. 
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(2) 


the brackets indicating time average. Thus, given the VAF the growth 
the second moment the diffusate can described. less information 
needed, the diffusivity along the axis can calculated from the relation 


This last based the agreement the detailed description 
(1) with the classical theory for times long compared with all relevant micro- 
processes (relaxation time, time between collisions, coherence time for 
external forces, cyclotron orbiting time, etc.). previous report (Bourret 
1960) diffusion equation was advanced which the VAF appears. This 
equation gives more complete description the diffusive process taking 
into account the VAF characterization the turbulence. not dealt 
with detail here, but Section III its use with the VAF shown illus- 


tratively. 


THE LANGEVIN EQUATION 


What wish calculate are the autocorrelation functions the velocity 
components representative diffusate particles. These autocorrelations are 
time correlations rather than space correlations. The latter are 
importance the theory turbulence and are calculated, after certain 
simplifying assumptions have been made, from the equations 
fluid motion from semiempirical theories. our present problem 
special simplifying assumptions need made once have decided upon 
the form the equations motion. Since, however, wish describe 
gaseous systems which the particles collide must abandon strictly 
deterministic equation and adopt, instead, the Langevin equation which 
describes not the exact behavior single particle but the behavior 
typical particle the ensemble. The choice this equation, however, imposes 
drastic limitations the kind information can obtain. The equation 


may written 


This says that the acceleration the average particle due applied 
force which may due collisions with other particles, the deflecting 
force the magnetic field and the retarding effect 
collisions with other particles with the medium which moves. The 
coefficient the retarding force, may approximately equated with the 
collision frequency. Its exact value, however, must calculated from con- 
sideration the details the collision process. (For detailed discussion 
this see Molmud (1959).) The force term, could specified explicitly, 
would make possible the calculation the trajectory the representative 
particle. order study the diffusion particles would necessary 


(3) = hm Riu lt')dt’. | 
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average over the equations motion for particles differing their initial 
conditions motion. that case would need know the particle density 
phase space order obtain the averages. What shall do, however, 
employ the Boltzmann equation for the distribution phase space 
basis for the derivation the Langevin equation. The distribution has thus 
been taken into account the formation the equation with which 
begin. shall consider, too, that the Langevin equation has been derived 
from the relaxation approximation the collision term 
the Boltzmann equation. The derivation from the Brownian motion approxi- 
mation stipulates conditions which cannot satisfied the cases interest 
the present inquiry (see, for example, Wang and Uhlenbeck 1956). the 
absence explicit force term must deal with the Langevin equation 
statistical fashion. This accomplished characterizing the force 
its power spectrum or, equivalently, its autocorrelation function. 
then possible determine the power spectrum the particle motions, and 
the conditions” the particle motion not appear. The only constant 
integration which requires specification this treatment the mean kinetic 
energy the particles. This equivalent specifying the temperature 
the diffusate gas, and assume here that the gas approximate thermal 
equilibrium all times. see, then, that the Langevin equation itself 
never solved when the force statistically characterized; only the velocity 
autocorrelations can deduced from it. possible, therefore, construct 
from the Langevin equation differential equation for the velocity autocorre- 
lations themselves, although generally proves easier construct integral 
expression for these autocorrelations directly rather than deal with the 
differential equations which they satisfy. 


VELOCITY AUTOCORRELATION FUNCTIONS 


First write the three dimensional Langevin equation for charged particles 
the presence magnetic field. index notation this reads 


010 

000 


The summation convention used here and the following. The term due 
the magnetic field has been written terms the matrix 

For simplicity the magnetic field supposed directed along the Since 
all the quantities used are tensorial character, the effect different choice 
magnetic-field direction can corrected for simply transforming the 
tensor quantities through appropriate rotation. The truncated Fourier 
transform the force components defined 
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analogous definition for the velocity components will used. 
assume for the moment that the force and velocity components vanish outside 
the interval —7/2 follows from equation (5) that 


where the mobility matrix, given 


(iw+g)[wo+ (iw+g) ] 0, 0, (iwt+g)° 


The velocity autocorrelation functions which are interested are defined 


1 +T7/2 
-T/2 
The set functions defined 


and called the cross-spectral density or, simply, power spectrum, related 
the velocity autocorrelations according the Wiener-Khinchine theorem 


Analogous equation (11), the power spectrum the force functions 
defined 


From equation (8) together with the definitions (11) and (13), follows 
that 


where the symbol has been introduced and defined 
will have occasion use the contracted expression 
which takes the form 


{ 
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where the auxiliary definitions 


2 


have been used. The velocity power spectra are given equation (14). The 
corresponding autocorrelation functions can expressed terms the 
quantities 


+o 
follows: 


Equations (14) and (20) are usually the best means obtaining the power 
spectra velocity autocorrelation functions the moving particles. 


ISOTROPIC WHITE NOISE FORCING 


The kind force usually treated connection with the Langevin equation 
thermal noise, and this may described isotropic white noise, i.e., the 
three directional components are completely uncorrelated with one another 
and each component has autocorrelation function representable delta 
function. terms spectral densities this means that there cross- 
spectral density between the components and that the spectral density 
single component white. This spectral density function represent 


(21) 


The constant has very immediate physical significance. Its meaning 
becomes clear from the normalization the velocity autocorrelation functions 
resulting when equation (21) employed equation (14). Because the 
definition (10) and the assumption thermal equilibrium the particle gas 
the velocity autocorrelations must normalized according 


The expression here refers the Boltzmann constant times the absolute 
temperature. This last should distinguished from the same symbol used 
previously different sense. The former usage will not repeated what 
follows and confusion should arise. This results the identification 
simple function the gas temperature through the relation 


(23) 


The autocorrelations the velocity components themselves prove be* 


*It may noted that this VAF corresponds the case given 
Bourret (1959). 
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The relation between the autocorrelation functions and the diffusion the 
particles has been described elsewhere (Bourret 1960). there proposed 
that the diffusion equation should be, the absence sources singularities, 


The two functions and are the polynomial numerator and denominator 
the Laplace transform the functions (24). 


wish apply equation (25) the description the diffusion charged 
particles perpendicular magnetic field first integrate over the co-ordi- 
nate (parallel the magnetic field) and introduce the definition 


The functions and are, two dimensions, 


(28) x(s) 
and the resulting diffusion equation 


The high order terms this equation describe transient effects which are 
little interest. any case they are usually small magnitude compared with 
the terms second order less. Neglecting terms order higher than the 
second equation (30) may written the form 


with 


This equation the form the Goldstein equation for turbulent diffusion 
and has been discussed the present author elsewhere. The diffusivity 
constant agrees with the standard results (see, for example, Allis 1956). 


CASE ISOTROPIC MARKOVIAN FORCING 


the forces are slightly more coherent than the previous case, and 
are, particular, Markovian, the particle velocity autocorrelations can 
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found from equations (12) and (14). shall assume that the applied forces 
are Markovian and isotropic and possess correlation functions 


and spectral functions 

(34) 


The resulting particle velocity autocorrelations prove 


where the notation 


(37) 


has been used for brevity. the limiting case which (36) assumes 
the form 


These autocorrelation functions are closely approximated those given 
equation (24) all cases likely physical interest—such the 
ionosphere. highly coherent radiation were responsible for the forcing, i.e., 
were not very small compared this approximation would longer 
valid. This possibility for laboratory situations. The cross-correlation 
terms have not been calculated here because they cannot contribute the 
diffusion according the proposed equation (25). This because the anti- 
symmetrical elements the Gag matrix (equation (17)) cancel the summa- 
tion the diffusion equation. The spectral functions corresponding expres- 
sions (35) and (36) are readily shown 


VI. THE CASE FORCING COLLISION IMPULSES 


the charged diffusate particles are mixed with much larger number 
neutral particles, such may occur slightly ionized gas, the driving 
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force will consist principally random sequence impulses due collisions 
with members the neutral species. The force correlation matrix defined 


where the brackets indicate ensemble average and the wavy line means 
average over all incident directions the collision. This last equivalent 
averaging the matrix over the rotation group three dimensions. The 
result this averaging simply eliminate the off-diagonal elements, i.e., 


(42) = 46, Fi(t’) +t) + Fo(t') Fo(t' +t) + F3(t’) F(t’ +0). 


Since the collisions are assumed isotropic each the bracketed terms will 
have the same average value and may write 


where typical force component. This force component consists 
random series positive and negative pulses functional form will 
assumed that the positive and negative pulses occur with average fre- 
quencies and respectively. The autocorrelation function appearing 
(43) can evaluated means generalization Campbell’s theorem 
(see Appendix) for random series positive and negative pulses which 
states 


2 
For the individual pulse assume the simple form 


(45) 
which the effective duration the collision and the total impulse. 
Substituting this pulse shape into equation (44) gives for the force autocorre- 
lation matrix (43) the expression 


will noticed that the second term equation (44) vanishes because the 


pulse series assumed have equal numbers positive and negative pulses. 
From (46) readily find the spectral matrix the forces 


2 —2rolw 
(47) ,;() = bu (Ap) 
From this the spectral matrix for the velocities found from equation (14) 


i 
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The classical diffusivities are given equation (3) and are given 


From equation (17) this becomes 


(50) 
equation (20) the VAF are determined the integral 


For the case with magnetic field this yields 


where the following notation has been used 


(53) I(t) 
(54) 


and where the normalizing constant has been determined the relation 


VII. CONCLUSION 


was mentioned the Introduction that velocity correlations space 
would not taken into account here. This omission makes impossible 
apply the results obtained plasmas undergoing collective oscillations since 
the nature these oscillations imply spatially coherent motion 
the charged particles. The discussion given here, moreover, does not include 
description oscillating particles subject random forcing, collision 
damping, and magnetic field deflections. Just how such situation might 
legitimately described Langevin equation not yet entirely clear. The 
effect the oscillation the VAF may prove very similar that 
the magnetic field, i.e., may simply add periodic factor the result. This, 
well the problem spatial VAF, will taken the author else- 
where. 
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APPENDIX 
wish construct the autocorrelation function series composed 
randomly occurring identical positive and negative pulses. For this purpose 
must construct slight generalization Campbell’s theorem, which, 
its ordinary form, applicable only series single kind pulse. The 
pulses will have the functional form f(¢), and subscript will indicate the 
algebraic sign the pulse, thus 


The entire series positive and negative pulses will designated F(t), 
and will written the sum series positive pulses plus series 
negative pulses follows: 


The and the are the random times occurrence the positive and 

negative pulses, respectively. The ensemble average used forming 

the autocorrelation function will indicated brackets. this notation 

the autocorrelation function for the series may written out sum 

four terms 


assume that the average numbers positive and negative pulses per unit 
time are and respectively. The form Campbell’s theorem given, for 
example, Rice (1954) allows write down the autocorrelation functions 
for the first two terms the preceding equation follows: 


The third term write out 


Carrying out the averaging over all possible values the and and 
indicating the average total number pulses per interval and 


Tc 
*This assumes, course, that the ensemble average may replaced time average. 
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and finally 
2 
since 
(66) 
T3020 T-+00 


The last term the complete expression for the autocorrelation the series 
gives the same result the above. Combining these formulas with the results 
given Campbell’s theorem for the first two terms equation (59) gives 
finally the desired generalization 
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NOTES 


THE SUDDEN INCREASE COSMIC RAY INTENSITY MAY 1960* 


sudden increase cosmic ray intensity occurred May approxi- 
mately 10:30 U.T. Previous events this type are rare and have been inter- 
preted burst particles coming from the sun the immediate neighbor- 
hood with energies sometimes few decades billion electron volts. The 
last event this particular type and the largest yet recorded occurred 
February 23, 1956, though unusual event July 17, 1959, had some 
similar features. The sudden increases the type that occurred May 
are believed always associated with intense solar flare but 
means occur with every intense flare. 

Four stations Canada are operated the Division Pure Physics 
the National Research Council. Each station equipped with cubical tele- 
scope (24 inches side) and standard neutron monitor. The stations are 
listed follows. 


Geographic Geomagnetic 
Name Lat. Long. (meters) Lat. Long. 
Sulphur Mountain 51.2°N. 115.6°W. 52.9° 
Ottawa 101 56.8° 351.4° 


The records from the four stations have been measured and the intensity 
variation shown Fig. The top four curves represent the changes 
the intensity the nucleon component measured standard neutron 
monitor each the above four stations. The 100% base represents the 
mean between from 0900 1030 G.M.T. each case. With our method 
recording, the pips can read off the recorder any time interval, though 
there increasing reading error for shorter intervals. The shortest interval 
used was minute. The reading error may represented the number 
counts +0.05 minutes though fact this error probably not normal 
statistical error and depends the reading technique used. This error should 
combined with the standard error for the Poisson distribution when the 
intervals are short but for intervals minutes more can neglected. 
The counting rates the 100% level are shown Table 

For the main curves Fig. the counting interval was varied from 
minutes depending the rate change. The inserted curves for Churchill, 


*Issued N.R.C. No. 5790. 
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Sulphur Mountain, and Ottawa show independent set measurements 
using minute intervals with expanded time scale. These are included 
show the shape the peak intensity. This well rounded each case, there 
being evidence appreciable higher counting rate even for few 
seconds. barometer corrections have been made for any these data but 
the barometric pressures were not changing rapidly enough affect the 
results appreciably. The lower curves represent the meson intensity change 
the four stations measured with standard cubical telescopes. The times 


TABLE 
Counting rate the 100% level counts per hour divided the scaling 
factor 

Nucleon component Meson component 

Scaling Scaling 
Station Pips/hour factor Pips/hour factor 
Resolute 220.7 51.2 1024 
Churchill 265.9 61.8 1024 
Sulphur Mountain 512 80.9 1024 
Ottawa 262.3 57.1 1024 


onset and times the peak are listed Table II. The times onset are 
estimated the point where straight line through the rapid increasing 
part the curve crosses the 100% axis and for the nucleon component are 
probably accurate better minute. The times maximum are 
estimated using certain amount smoothing the minute values but 
cannot error more than about minutes. the case the meson 
component the accuracy lower because the much lower amplitude 


the event. 


TABLE 
Time onset and time and height the maximum intensity 
Station onset maximum maximum 


Nucleon component 


Resolute 1032 
Churchill 1030 1037 390 
Sulphur Mountain 1029 1042 420 
Ottawa 1029 1040 298 
Meson component 
Resolute 100 
Churchill 1028 1036 113 
Sulphur Mountain 1031 1038 116 
Ottawa 1028 1035 108 


The striking feature the results are: 
(1) the low value the increase Resolute relation the other stations. 
The increase either lacking insignificant the case the meson com- 


ponent. 


ns. 
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(2) the time delay reaching maximum Resolute. This quite signifi- 
cant, the maximum being least minutes later than other stations. 

(3) the high intensity Churchill relative Ottawa and Resolute, all 
which are low altitude. 

(4) the rise Churchill seems steeper and the maximum slightly 
earlier than Sulphur Mountain and Ottawa. 

(5) the high intensity Sulphur Mountain, which may least partly 
due its altitude. Sulphur Mountain and Ottawa have about the same 
geomagnetic latitude. 

(6) the relative heights the nucleon and meson components are roughly 
consistent the different stations. 

Before drawing definite conclusions about this increase, the data from other 
stations should studied but certain preliminary qualitative conclusions can 
presented. 

consideration calculations Liist (1957), none these stations were 
impact zones the time this increase, though Sulphur Mountain may 
have been not very far from one. The declination the sun was +16°. 

The flare reported to.be associated with this increase was very short 
duration, and, according the preliminary listing, started well before the 
cosmic ray increase, namely, 1020 (Preliminary Report Solar Activity 
from the High Altitude Observatory, Boulder, Colorado, 453 May 
1960) and was over 1025. The times onset Table indicate that 
there was delay between and minutes the arrival particles. This 
delay the same order that found the February event, 1956, for 
non-impact zone stations (Liist and Simpson 1957). 

The small intensity Resolute suggests that the isotropy suggested 
Meyer, Parker, and Simpson (1956) was not achieved during the major part 
the event. superposition the part the increase after about 1110 
shows close resemblance between the curves for Resolute, Churchill, and 
Ottawa, indicating that isotropy was probably reached about then 
minutes after onset. 

The data from which the curves Fig. are drawn have been distributed 
World Data Centers and the I.G.Y. cosmic ray stations with which data 
exchange has been arranged for some time. Spare copies are available. 


The author would like express appreciation the prompt forwarding 
data Dr. Terentiuk from Sulphur Mountain, Mr. Gustafson 
from Churchill, and Mr. Bridges from Resolute. 


Ltst, 1957. Phys. Rev. 105, 1827. 
and 1957. Phys. Rev. 108, 1563. 
P., PARKER, N., and 1956. Phys. Rev. 104, 768. 


RECEIVED JUNE 21, 1960. 
PURE 
NATIONAL RESEARCH COUNCIL, 
CANADA. 


its 


LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should submitted not later than six weeks 
previous the first day the month proof will sent the authors. 


Frequency Measurement Standard Frequency 


Measurements are made Ottawa, Canada, using N.R.C. caesium-beam frequency 
reference standard (with assumed frequency 192 631 Frequency 
deviations from nominal are quoted parts per negative sign indicates that the 
frequency below nominal. 


GBR, 
Date, 
June MSF, 5-hour 24-hour WWVB 
—152 —150 N.M. 
Midmonthly mean — 156 —150 -—150 


Midmonthly mean 
WWV 


Note: N.M. no measurement. 
*Time of observations: 01.30 to 03.30 U.T.; 10.00 to 13.00 U.T. 
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Ottawa, CANADA. 


1Issued as N.R.C. No. 5783. 

*Cf. Kalra, S. N. 1959. Can. J. Phys. 37, 1328. . 

4Values published previous to this issue were subject to an additional uncertainty due to the method of assessment 
of the magnetic field correction used in the resonator. 
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text. Titles should always given but should brief; column headings should brief and 
descriptive matter the tables confined minimum. Vertical rules should not used. 
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General.—All figures (including each figure the plates) should numbered con- 
secutively from up, arabic numerals, and each figure referred the text. The 
name, title the paper, and figure number should written the lower left corner the 
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should not more than times the size the desired reproduction. Whenever possible 
two more drawings should grouped reduce the number cuts required. such groups 
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made for the captions. The original drawings and one set clear copies (e.g. small 
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Photographs.—Prints should made glossy paper, with strong contrasts. They 
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duplicate; they are reproduced groups one set should mounted, the duplicate 
set unmounted. 
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